For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


Gx pais 
UNIOTRSTTATS 


The University of Alberta 
Printing Department 
Edmonton, Alberta 


Digitized by the Internet Archive 
in 2024 with funding from 
University of Alberta Library 


https://archive.org/details/Vossoughi1973 


THE UNIVERSITY OF ALBERTA 


NON-NEWTONIAN FLOW DISTRIBUTION IN A POROUS BED 


by 


"aaa VOSSOUGHI 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH IN 
PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF 


MASTER OF SCIENCE 


in 


CHEMICAL ENGINEERING 


DEPARTMENT OF CHEMICAL ENGINEERING 


EDMONTON, ALBERTA 


FALL, 1973 


Nel Aahjove! 


whom I Love most 


iv 


ABSTRACT 


Fluid flow of both Newtonian and viscoelastic 
fluid was studied in a model porous medium composed of a 
matrix of small glass cylinders. Local velocity profiles and 
flow distribution were determined by a streak photomicro- 
graphy technique. Measurements were obtained for average 
interstitial velocities of the bed ranging from V/e = 0.0568 
cm/sec to 0.568 cm/sec. The major effect focused on was the 
variation in flow between pores owing to small inhomogeneities 
in the bed. 

Compared to the Newtonian fluid, the effect of 
viscoelastic fluid was found to be a suppression of the 
velocity and also flow distribution variations with pore 
opening. The effect of elasticity on flow distribution was 
observed to occur before its influence in pressure drop 
measurement. 

Qualitative visual observations also revealed that 
the number of dead zones, or nearly dead zones and recircul- 
ation regions for viscoelastic fluid were reduced relative 


to the Newtonian case. 
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Re 


At 


Total volumetric flow rate, em /4ec. 
volumetric bes rate passing through each 
opening, om?/sec. 

One half the gap between two parallel 
plates, cm. 

Reynolds number defined by Equation (V-9), 
dimensionless. 

Reynolds number defined by Equation (B-4), 
dimensionless. 

Gas constant, gr om? /sec* g-mole oK 
Hydraulic radius, cm. 

Radius of the platen, cm. 

Revolutions per minute. 

length of the streaks, cm. 

Length of the streaks, in. 

Sensitivity of the torque transducer, 
muUcron/ Voltz 

Sensitivity of the normal force transducer, 
Gu/N7o le. 

Minimum free flow area, cm”. 

longitudinal piteh, sem. 

Specific surface, ae 

Transverse pitch, cm. 

Wetted surface per volume of bed, oe 
Duration of one revolution, second/rev 


Movement of torsion head transducer, volt. 
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Absolute temperature, oK 


Time of period of opening or closing the 


light, sec. 
Torque, dyne cm. 
Interstitial velocity, cm/sec. 


Approach velocity, cm/sec. 


Average interstitial velocity from center 


line openings, cm/sec. 

Average value of the maximum velocities 
center line openings of two consecutive 
cm/sec. 

velocity in the z-direction, cm/sec 
Superficial velocity, cm/sec. 

Average velocity of the flow through 
parallel plates, cm/sec. 

Volume occupied by cylinders, com>. 
Output of transducer, volt. 

Total volume, com>. 

Void volume, com>. 

Volume of the water, com>. 

Mass flow rate, gr/sec. 

x-coordinate, cm. 

X-coordinate, in. 

y-coordinate, cm. 


z-coordinate, cm. 
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Greek Symbols 


« Angle of Panapraiotess 

B Angular rotation of the platen, rad./sec. 
V Delta operator. 

€ Volumetric porosity, dimensionless. 

E. Surface porosity, dimensionless. 

ik Average nominal shear rate, sec +. 

2 Mobility, cm? sec/gr 

u Viscosity, poise. 

Wo Limiting viscosity at low shear rate, poise. 
Ue Viscosity of the solvent, poise. 

ie Density, gr/com>. 

L Shear stress, dyne/cm 

0 3 Limit of the relaxation time at low shear 


rate, sec. 


O- Relaxation time of fluid, sec. 
opr Duration of a process, sec. 
0 Shear rate, gectek 
Subscripts 
B indicate an ideal bundle of tubes. 
e elastic 
G Glycerol 
fe) On 
S Separan 
WwW water 
W wall 
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modified for the porous bed. 
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CHAPTER I 
INTRODUCTION AND BACKGROUND 


Lyk General 

Increasing demands for petroleum supply and the 
limited number of the oil reservoirs make it inevitable not 
only to use the produced oil efficiently but also to increase 
the ultimate recovery of reservoirs. Investigation of the 
many fields per eeed (1) that more than 50 per cent and in 
the case of heavy oil reservoir, even 80 per cent of the 
Original oil in the reservoir remained unrecovered after 
primary recovery; i.e., when the production rate became 
so low as to make economical operation impracticable. 

Oil production requires energy to expell the oil 
trapped in the interstices of the porous rock into the pro- 
duction well. This energy may be available in the reservoir 
initially in different forms. The major forms of natural 
energy in the oil reservoirs are: 

1. Compression energy of the oil or water within the 
field or contiguous to the oil-bearing rock. 

2. Compression energy of the gas dissolved in the oil 
or free-gas zones overlying the oil layer. 

3. Gravitationali energy, due jto? the difference..in hydro- 


static level of producing well and upper part of the 
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reservoir. 

If no natural energy is available in the reservoir 
or if the existing energy is not sufficient, the external 
energy should be applied to exceed the economical rate of 
production. Even in the case of high pressure reservoirs, 
the initial content of energy would be depleted after a long 
period of production and, therefore, external energy would 
be required. 

Fluid injection as an external energy supply to 
the reservoir may be applied for two different purposes; 
(a) during the period of the operation of the well to com- 
pensate for the pressure decline created due to the oil 
withdrawals from the space voidage. The fluid injection 
will help to keep the rate of production almost constant 
and above the limit of economical rate, (b) after the pro- 
duction well reaches a state of substantially complete 
depletion of its initial content of energy. In reservoir 
engineering terminology the former is called "Pressure 
Maintenance" or "Pressure Regulation" and the latter 
"Secondary Recovery." In secondary recovery operations the 
pressure required for the fluid injection is often lower 
than that which would be required for injection under 
primary pressure maintenance. 

Fluid injection as a method of oil recovery is 
indeed more than one century old. The first accidental or 
perhaps intentional water flooding occurred in 1865 in the 


Pithole City of Pennsylvania (2) Since then a great variety 
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of fluids such as gas, liquids, even slurries, and also 
different additives to the fluids to achieve better per- 
formance have been tried. 

Of the different methods of the secondary oil 
recovery, water flooding is the most advanced method. 
Gas-injection operations, usually termed "gas repres- 
suring," are less common than water flooding. 
Secondary recovery by water flooding has been carried on, 
aS a common practice, in regular patterns. Historically, 
the first type of regular network was the alternating-line- 
drive pattern(3). Next followed the five-spot pattern which 
is most commonly employed. The mechanisms of the fluid 
displacement for various patterns of injection and produc- 
tion wells have been extensively investigated (4-12). The 
injected water pushes the oil in a somehow piston-like 
behavior from the vicinity Of the injection wells and ahead 
of the advancing water, a zone of high oil saturation is 
formed which moves toward the producing well. The clean 
oil is produced until the oil-water bank reaches the producing 
well at which point breakthrough occurs. After breakthrough, 
increasing amounts of water are produced with the oil. 

Development of the oil bank and its stability 
governs the success of the operation. In some projects 
the high oil saturation zone does not develop at all and 
breakthrough takes place from the start. The rate of the 


production of water may reach a limit such that the production 
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of oil becomes too low to support continued operation. 

It 18 intuitively understood that if the conditions 
at the start of the flood are such that the water can move 
more readily than oil, the displacement front will spread 
into the oil and early breakthrough occurs. Water intrusion 


into the frontal oil bank is termed "fingering." 


i pe Mobility Ratio: Concept 

Whether or not a stable oil bank develops depends 
largely on the relative permeability characteristics of the 
formation, on the viscosity of the displacing and displaced 
fluids, on the saturations of water and oil at the time the 
water injection is started and also on the inhomogeneities 
of the rock itself. The first two conditions can be regrouped 
in mobility ratio (13) which is defined as the ratio between 


the mobilities of the displacing and displaced fluids, where 


(LL) 


m represents mobility ratio and ee and 2 are mobilities of 
water and oil respectively. 


Definition of the mobility is better understood by 
studying Darcy's law. The theory of flow through porous 
media originated by Darcy's experiment in 1856 (14a). Darcy 
found that the flow rate is proportional to the pressure 


drop along the bed, or in terms of superficial velocity, 


i 


rox J 4 hat ee uti bs £ fesiZ ' 
ae a 
vert eo Ysew eds 4 re’ on io a Yo a 
eM (fiw trond scuseianal qekb at vito nest ott 
nore Ye7 eW 6S 7LD90 onsite pad etn be 
“pri stepak2" home ek Aurea Lio nant wit O8f 
— _ 
e = 
. +qe0deD oka +L Lido 
a = 
i» aeqgoleveb Aned 160 elds In ‘en 146 seat 
) Xo 122 tevio yt dssigtsu Suisse tax ont ee 
, 
hex rib bas paboelquis ett to ytieooety ama ae of 
a 
; wer x 14 Ss if © f Pee.) y y mw ¢2 e "Of3 iD T4568; ena ie 
; 2s 
Asis ti ott ao obis bas bedtbse |r ioksoatnk 
) 291 od ago enotdibroo ows seuci edt .TIeedt tao 
f ose 44 boanttuh 2f duidw (£1) cigem veto 
“4 ebiol? beosl¢elb bas peices iqeth ont To eoktt tide 
r a” 4 j p 
f Ww | y" 
f i ' — 2 
oo : 
« ; ~~. 
i 4 LJ de "isis “J ; . ‘ byris Orvgs2 yot Edom eatnses Awl. @ 
t) ; q 7 ied ts 
7 "| ‘ A Pi 
a é fevi 4s 5 8 fio Gis” (936 
a 
enabat wetted ef vai (rdom‘sdiy Yo Aotsbarts y 
; : ae = 
avoszeq dovetidd wok? 20 ¢rosds ent ..wal a3 160 pare ‘bist 
| ; Stes wey 
yous . (eb) a@@l al taonc-: 19QXS a 'yousd vd be we a] 
’ x, 
sxyeeerta ofy of Lecotidsometq et eae? wold ent 


wiitocoliev Isiotiaeawe to amiat stk x10 rs dt pene 


a 


AP, 


V=A L (1-2) 


where V is superficial velocity, AP, pressure drop, L, 
length of the porous medium and i, constant which by defini- 
tion is the mobility of the fluid. Later, in 1930, Nutting 
elucidated (14a) the physical significance of \X and stated 
ial at & ie 


A = ran aay (I-3) 


where u is the viscosity of the fluid and K the specific 
permeability of the porous medium. Substituting (I-3) into 
(I=2), ‘one obtains, 


eh CO ae in 
| aw a 3 (I-4) 


Equation (I-4) is known as Darcy's law which holds for the 
laminar flow of Newtonian fluids through homogeneous porous 
media. 

Darcy's law was modified by Christopher, et al. (15) 
for the flow of power-law model fluids through porous media. 
Using a modified Blake-Kozeny equation, they obtained, 


Howe , 


Gi )) 
where n is power-law parameter and H, viscosity level para- 
meter, a factor which accounts for the additional dependence 
of V on K and porosity due to non-Newtonian behavior. 

The first theoretical attempt to show the role of 
the mobility ratio in fluid displacement was carried on by 


Leverett, who published his fractional flow formula in 1940. 


On neglecting capillary pressure gradient and gravity effects 
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his equation takes on the following simple form (l6a) 


fk OS GD) Ac 
WwW K (I-6) 
Léereanontw 
K 
w Oo 


or, uSing the mobilities of two fluids, 


——— ai 
1 + e) 
— 
WwW 
where, 
ae g im. 
ty 4 we : (I-8) 
W O 


q, and ae are volumetric flow rates of water and oil respec- 
tively. In an attempt to calculate the oil recovery, Welge 
(17), Aronofsky (4,5) and other investigators (8,12,18) also 
recognized that through decreased mobility ratio an improve- 
ment in displacement efficiency could be obtained. Mungan (19) 
amongst others showed qualitatively the conditions of stability 
at a displacement front in a simplified model. He considered 
the incompressible horizontal linear displacement of oil by 
water in a piston-like manner and neglected capillary effects. 
In a system of length L the displacement front is at position 


x where x < L. Writing Darcy's law for the positions filled 


with water and oil, 
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where q is volumetric flow rate, A, cross-sectional area, 


2, velocity of the displacement front, Pir Pu ger 


x a Gist 
at inlet, displacement front and outlet respectively. Elim- 


pressure 


inating 2b between Equations (I-8, 9), 


De eee ‘ 1 ° 


O W LA, ne x(A XQ) 


The Equation (I-10) shows, qualitatively, the role of the 
mobility ratio in stability of oil sweeping. When LV ee A 
i.e. he? Arg! for any increment increase of x, say dx, which 
physically is indicated by forming a tiny finger or kink of 
water into the oil-filled segment the velocity of displace- 
ment front, i.e. gqg/A, will decrease. Therefore, the finger 

or kink of water can not proceed into the oil-filled segment. 
That is the tendency for early breakthrough will be suppressed. 
For the case when haf wg <te l.e. eau any incremental increase 
of x will increase the velocity of the displacement front 
which leads eventually to complete breakthrough. 

Slobod, et al (10, 11) studied fingering in miscible- 
phase displacement by an X-ray shadowgraph technique. Their 
photographs showed that at high flow rate (i.e. 29 ft/day) 

a large number of very thin fingers was developed; whereas, 
the slow flowrate (i.e. 1.6 ft/day) showed a gradual blending 


of fluids with no fingering. Both had an unfavorable mobility 


Fatico, of @. )Photographs taken for: mobility ratio of unity 
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revealed no fingering. One interesting effect was the 
complete elimination of fingering when they used a continuous 
Gradation, in Viscosity of flooding fluid. It as worthwhile 
to note that the overall viscosity ratio had the unfavorable 
value of eae = 3. This can well be in support of Mungan's 
modification in polymer slug process (20). 

Since the mobility of the oil is a fixed parameter-- 
depending on the kind of reservoir and oil--the only possi- 
bility to maintain stability in water flooding is to reduce 
the water mobility. 

An improved mobility ratio would also increase 
recovery efficiency by correcting permeability distribution 
problems as discussed by Stiles (21), Dyskstra, et al. (22), 
and Pye (23). 

Mobility reduction may be achieved by one or a 
combination of the following two ways: 

1. Permeability reduction of the porous media. 

2. Increasing the viscosity of the injected water. 
Materials such as Glycerol, sugar, starch, etc. have been 
used to increase the water viscosity. However, due to the 
large amount of the materials to be used, the method is not 
economically feasible. The concept of using natural high 
molecular weight polymers to increase the water viscosity is 
not new. The first patent appeared in 1944 by Detting (Shell 
Development Co.). The reader is referred to 27 different 


patents mentioned in reference (1). These patents were 
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granted during 1944 to 1961. 


i. 23. Polymer Flooding 

The use of synthetic polymers in the injected water 
of the water flood, as compared to the former mentioned addi- 
tives, is a fairly new technique. Polyacrylamide was first 
introduced in 1964 by Pye (23) and Sandiford (1). Their 
experiments showed that the presence of this kind of polymer 
in dilute concentrations decreases the water mobility 5 to 
20 times more than would be expected from measurement of the 
solution viscosity. Since then it has become of increasing 
interest to study and investigate the flow behavior and 
mechanism of polymer flooding and recovery (20, 24-28). In 
addition to the laboratory experiments, the polymer water 
flood has also been successful in the majority of the small 
or even large scale field projects where it has been employed 
(29-3 3.e 

In polymer water flooding a water-soluble high 
molecular weight polymer which is partially hydrolyzed is 
desirable. The polymer solutions, at the concentrations 
generally used in polymer floods, are non-Newtonian and show 
pseudoplastic behavior in viscosity measurement and elastic 
effects in normal force measurement. The pseudoplasticity 
is suppressed greatly in the presence of the strong electro- 
lytes (20, 24, 34). 

The main preference of the polymer water flooding 


over the conventional method is that a greater fraction of 
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10 
the reservoir volume may be swept at water breakthrough. In 
other words, the sweep efficiency increases by introducing a 
small amount of polymer into the injected water. In addition 
to the increased macroscopic sweep efficiency, improved micro- 
scopic displacement at a given amount of fluid injected has 
also been observed (1), i.e. the residual oil saturation would 
be less in swept areas compared with that of conventional 
flooding. The surprising ability of the polymer to lower the 
mobility of flooding water may be attributed to the adsorption 
of the polymer molecules on the pore surface, to the plugging 
of the very fine pores by polymers, to the viscoelastic 
effect of the polymer solution, or to a combination of all 
of these. 

Lee, et al. (27) studied the behavior of the polymer 
solution in oil recovery by looking at pseudoplasticity effect. 
They concluded that provided the polymer concentration in the 
flood water is such that the average mobility ratio is favor- 
able or unfavorable, the oil-water contact may exhibit a 
stable or unstable situation during the course of the flood. 
Their visual study on a Hele-Shaw model verified their 
qualitative prediction of the fingering based on pseudo- 
plasticity effect. However, there is considerable evidence 
to doubt that the power-law model is sufficient to predict 
the flow behavior of polymer solution in real porous media 
under all conditions of interest. 


Burcik (35-38) attributed the mobility control of 
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the polymer solution not only to the viscosity improvement 
but also to the reduction of ‘aie effective permeability 

due to adsorption of polymer molecules on the surface of 
pores. He observed dilatant behavior of polymer solution 
flow through a disc of Berea Sandstone which, presumably, 

was a result of the bound polymer molecules within the pore 
structure reducing the effective permeability. To distinguish 
this behavior from true dilatancy, the author suggested that 
the phenomenon observed be designated as pseudo-dilatancy. 
Plugging effects, as well as adsorption have also been noticed 
bya other” investigatorse(205t 24) 28 29,931). tn another 
paper, Burcik and Walrond (39) proposed that small particles 
of microgel might be primarily responsible for the perme- 
ability reduction; however, Lynch and MacWilliams (40) in 
response to Burcik's claim showed that gel particles need 

not be present to have polymer units of the order of one 
micron and single polymer molecules of that size can be 
reasonably expected in the ranges of 3 million molecular 
weight. Patton, etal. (41) revealed that the polymers with 
less adsorption affinity on the pore surface are more effec- 
tivelatithe floods front’ ands alsovon the o11 recovery.. They 
found that the quantity of adsorption varied widely from one 
polymer to another. For one specific kind of polymer the 
adsorption varied from one type of mineral surface to another 
and increased with salt concentration (28). Reduction in 
relative permeability observed by Patton,et al. (41) was 


attributed largely to the tight connate water-saturated 
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pores. These colloidal particles could physically block 
the smaller capillaries and inhibit the displacement of the 
connate water. 

The relative permeability becomes insensitive to 
the concentration of the polymer solution above a minimum 
value of concentration. The dependency of the relative 
permeability on the polymer concentration has been discussed 
by Jennings, et al. (34). 

At this stage it is felt fruitful to define two 
factors namely "Resistance Factor" and "Residual Resistance 
Factor." Both are descriptive of polymer behavior but in 
different situations. Resistance factor, as defined by 
Pye (23), is the ratio of the brine mobility to the polymer 
solution mobility. Resistance factor is indeed a measure 
of adsorption together with entrapment of the large mole- 
cules in the small size pores and also elasticity which 
mostly has been ignored. Jennings, etal. (34) showed that 
resistance factor decreases as K/e increases; where K is 
permeability and e«, the porosity of bed. Residual Resistance 
Factor called "Permeability Reduction" by Smith (28) is 
defined as the ratio of the brine mobility in a rock, not 
previously contacted by polymer, to the brine mobility in 
the same rock after being swept by polymer solution. Smith 
(28) pointed out that permeability reduction is due to the 
immobile polymer retained within the porous medium. It is 


greater in less permeable rocks and increases with the 
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13 
polymer molecular weight and also with the flow rate. 
According to Jewett, et al. (32), the oil recovery increases 
Wat both facters* which Ws "nconsistent* with Patton, et ail. 
(A1)% 

In flows which are non-steady from a Lagrangian 
point of view, the flow behavior of polymer solutions is 
often surprisingly different from that of Newtonian fluids 
and is very dependent upon exact details of velocity field; 
e.g. flow around submerged objects or flow in converging- 
diverging channels. Ultman, et al. (42), studied the flow 
of polymer solution around a submerged single cylinder and 
sphere using Polyacrylamide and high molecular weight 
Cellulose Gum as polymers. They introduced Weissenberg 
number and studied the effect of the magnitude of the 
Weissenberg number on the effect of viscoelasticity. 
Weissenberg number was defined (42) as, 

Ceall 


aes = DT CLs) 


where Ce is the limit of relaxation time at low shear rate, 
U, the approach velocity, and D, diameter of the cylinder 

or, sphere. They noticed a significant difference. in the 
shape of the streamlines. For instance the streamlines of 
the flow of viscoelastic fluid across a submerged cylinder 
responded to the presence of the cylinder earlier than those 
of the Newtonian fluids and returned to the undisturbed 
motion sooner. The position of the separation point at the 


aft portion of the cylinder was noticed to be dependent upon 
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the magnitude of the Weissenberg number. 

The effects of elasticity of polymer solutions in 
porous media have been completely or partially ignored in 
a great number of papers published in this field. Some 
authors such as Jennings, et al. (34), Smith (38), and 
Harvey (43) believed that the effect of elasticity of the 
low polymer concentrations of polymer solutions in porous 
media is insignificant. Sadowski and Bird (44), to justify 
their capillary-hydraulic radius model based on three- 
parameter Ellis model which suffered from lack of the time- 
dependent elastic phenomena, reasoned that no significant 
elastic effect would be observed provided the fluid relaxation 
time is small with respect to the transit time. 

By the same reasoning one can conclude that the 
effect of elasticity could be significant if the geometry 
of the media is such that the transit time of the fluid flow 
through a contraction or expansion in a tortuous channel is 
comparable with the relaxation time of the fluid. This 
expectation has been theoretically shown and experimentally 
verified by Marshall and Metzner (45). They gave a correlation 
of the effect with Deborah number for Polyisobutylene and 
partially hydrolyzed Polyacrylamide solutions. Deborah 


number was defined as, 


N Bh eek coals 


where 0. is the relaxation time of the fluid, and ont the 
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1b) 
duration of a process which represents either the duration 
of a given deformation state or, equivalently, the reciprocal 
of the rate at which the deformation process changes (46, 47). 
In their analysis, Marshall and Metzner (45) used the convected 
Maxwell model to approximate the behavior of the viscoelastic 
fluid and thus could estimate the fluid relaxation time from 
steady shear flow measurements of the normal stresses. By 
assuming the porous media to be a series of converging 
sections in the flow direction Z, they showed that the 
appropriate measure of the process time in a packed bed of 


spheres is the reciprocal of the deformation or stretch rate 


du 
inthe Pilow “direction; ive: i) 1 


mately D,p/4, where Un is average interstitial velocity and 


which in turn is approxi- 


Dp the particle diameter. With this choice of Or! Equation 


(I-12) becomes, 


ou, 6 
Uepaae se eee (I-13) 


The quantity “m/D, may also be viewed as a measure of the 
deformation rate or so called elongation rate in the flow 
direction. In pure steady elongation flows, the Maxwell 
model they used predicts infinite stresses must occur when 


the N defined by (I-13), approaches 1/4. Intuitively 


Deb’ 
then, they suggested that the pressure drop required to 

pump a viscoelastic fluid through porous medium would rise 
to infinity as Nneb defined by (I-13) approaches a value of 


1/4. The experimental pressure drop measurements presented 


by Marshall and Metzner indicated that appreciable influences 
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of the fluid elasticity were felt when N as defined by 


Deb’ 
Equation (I-13), was approximately 0.05-0.06 and at the very 
low values of the Deborah number the effect of elasticity 

were negligible. It is worthwhile to mention that some 
investigators (15, 48) have failed to see the viscoelastic 
deviation in friction factor--Reynolds number plot even at 
higher values of Deborah number than that from Marshal and 
Metzner (45) or Sadowski (44). 

Elasticity has no effect on flow behavior in 
Lagrangian steady flows. Therefore, the increase in pressure 
drop due to elasticity can only be observed in an elongational 
flow experiment; i.e. a converging-diverging strata which is 
very much likely to appear in porous media. The failure of 
the capillary tube model of porous media to predict this 
phenomenon is, therefore, evident. 

Pye (23) seems to be the first one who observed the 
increase in resistance factor with advance rate which appeared 
to be noticeable only in tortuous passages not straight 
capillaries. Similar observations and statements have been 
brought up by Gogarty (24), Dauben and Menzie (25), Smith 
(28), Jennings ,~et.al.6(34) fandvHarrington,. etials (49). 1t 
is interesting to mention that Jones and Maddock (50) observed 
this behavior of viscoelastic fluid in porous media, but they 
considered it disadvantageous rather than advantageous for 
oil displacement. 

The effect of elasticity of the viscoelastic fluids, 


contrary to what Jennings, et al. (34) stated, need not, 


- bent? tab 26) 


d% 


yrov stilt de Hos 30, O- ~2Oy4 


ytiviseslo te. dontteul 
* be 


mon | >) ; 
all oe hes 
gnoe teats cotstaen ot @ zOwW 


idasleovatvy orld see og pe hae ovat 

Ly : - 7 * 

ib Mmeve FOL TAG tts im = Lomyese-3 LOdoeR HO, nti a 
T 3 Sieg 


bie L[edevsM mous? Jang apis. tecue slerods a. ao 


. g. att i 
. ‘i (oe, Hewobss * 


yoftvaedsed w > Ko toate on eset ok hanes 
a « 


GSIGDeaRia f t aseagonu0at odo ,eseieveta meal? be te, 


i | ° : : , an 7 “wy weap ——— 7 — 
hee | ' Tae So tad <9 . : ve ” & fi! fv er > 


to suulis3} off .esibem seoteg mE teenge od, ee 
sict Soibexq.ot sibemoevosoqg Jo [6 om edu coats 


.jnabive edo t08 otk sth ON 

‘ : 7 . — ; rhe 
it beyisado ofw eno sexi? off ed os |emope sy in eee, Ka 
Dal ssags fioisw estar eonevBe div wots ost sonedateon mh - 
iopisuse jon sepssesy svondaes at vilinio auntie r: a 
it a : oy 

iP a ny 
te bas enoljievresdo xacbathe ~ es b ae It 


cia) 7 


a 


iteme ,(¢8) etsaeM Bae mor anal «fb "s) vd tsbOD* Di: ae a 


Ji 6. (88) «fs #5 .nospnizias bose (26) »e ag 
DoevIecedas. (Ud) 120 DbhoM DTLS 4anou Seis aoishem od a — ag ae" 
' 
= r 4 4 Z _ - . ” * 4 77 E < 7 if _— _ abi 
yais stud ,sibem svortoqg at beuild Gs senleo soosty Fo‘ * 


7 a? 


aa a“ 


,86iull oitesloooriv edz oe ) bes ts 


.toa been besase (wey. ‘te’ é% 


how i 


i 


: - 


ery 


17 
necessarily, be expected to happen only at high flow 
rates. Jennings, et al. isus@imaea to include the relaxation 
time into the flow resistance correlation. They considered 
the total pressure drop, P, across a system to be the sum 
of the viscous, kinetic energy, and an elastic term, 
AP ot %. Oy 


AP a= Cy Hats VonteGs oruntt AP or (I-14) 


where Cy and Cc. are constants and 9, density. They made a 
further assumption that the elastic pressure drop be pro- 
portional to the normal stress difference. The final 


correlation appeared to be, 


AP On 
ee te Co ee thea Ul 0 ——5" (i= 15) 
m d 
S 
where, C, is constant and das the diameter of the constric- 
tion. They plotted Am V.S. u, for two different molecular 
m 


weight polymer solutions passing through a core specimen of 
rock and noticed a significant difference in the slope which 
was presumed to be due to the sensitivity of 9®f--not visco- 
sity--with molecular weight, and quantitative justification 
was avoided. 

To conclude, this entire section can be summarized 
as follows. Addition of polymers to water has been shown to 
be effective in enhancing behavior of water floods in porous 
media. The mechanism and therefore EoReEoliing factors, 
however, are poorly understood. In a number of experiments 


increases in pressure drop, much above what would be expected 
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18 
from a consideration of viscosity, have been observed. This 
effect has been variously attributed to pore blockage, adsorp- 
tion and viscoelasticity. In the case of viscoelasticity, 
theoretical considerations, although only semi quantitative, 
suggest significant changes can occur in flow distribution 
and in pressure drop. It is significant that the predictions 
and resultant dimensionless groups are very dependent on the 
velocity field of the fluid. As pointed out, in steady flow 
in a capillary tube, elasticity has nothing to do with either 
pressure drop or velocity distribution. 

Since experimental data and an analysis are available 
for flow around a single cylinder, it was decided appropriate 
to study the flow of a viscoelastic fluid in a porous bed 
composed of closely spaced cylinders. The scale of the 
experiments is such that any effects of blockage, etc. can 
be safely ignored and attention can be focused on any changes 
in flow field that might occur as a result of continuum pro-~ 


perties of the polymer. 
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CHAPTER t2L 
THEORY 


aE. al telintroduction 

Solving the equations of change to predict the 
velocity field and pressure drop with the exact boundary 
conditions of the present model would be extremely tedious-- 
if at all possible. However, as a first approximation, 
for slow flow the opening between two parallel cylinders-- 
especially when the gap between is small compared to the 
diameter--can be simulated by two parallel plates with an 
unknown equivalent length. The unknown equivalent length for 
the entire bed can be obtained experimentally by experiments 
with a Newtonian fluid. The procedure would be the same as 
for the capillary tube model commonly used in porous media. 

The capillary tube model, in spite of its lack of 
ability to portray time-dependency, has been successfully 
applied in porous media to predict the pressure drop or 
friction factor--Reynolds number relation of purely viscous 
flow (15, 51). The model is based on cylindrical Poiseuille 
flow using hydraulic radius of the medium instead of radius 
of the cylinder. Then, a correction factor due to the 
tortuousity of the bed is considered. A similar approach 
but with parallel-plates approximating the capillary pores 
will be followed for the present bed. It is worthwhile to 
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20 
mention that the pores of the present bed are actually 
rectangular channels with dimensions d and &, where d is 
the opening between two rods and &, the exposed length of 
rods. However, since £ >> d the channel can be fairly taken 
as a two infinite parallel plates. 

To achieve a more comprehensive feeling concerning 
the above approach, consider the plane perpendicular to the 
paper along the line A-A in Figure (II-l). The space between 
the glass rods along the line A-A are simulated as parallel 
plates and at the rear side of the rods, say along the line 
B-B, the plates are discontinued. Since the gaps are very 
small most of the pressure drop will result from the flow 
between adjacent cylinders in a given row rather than flow 
from one row to the next. 

The shape of the intersticial velocity profiles of 
the assumed model is the same as that of flow through parallel 
plates. The equations of the velocity profiles for the flow 
of Newtonian and power-law model fluids through parallel 
plates, which can be found elsewhere (52a, 53) and are 
presented in Appendix D, will not be discussed in this 


chapter. 


II. 2 Parallel-Plates Capillary Model--Pressure Drop 
Consider the flow of a fluid through two parallel 


plates in the z-direction as shown in Figure (II-2). One 


can write, 


qo = 27 fu, dy dx, ore 
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23 
where, qy is volumetric flow rate, uv velocity in the 
z-direction, %, the length of the channel in transverse 
direction and R, half the gap between the plates, i.e., 

Ryw= d/2- (II-2) 
Since flow is assumed rectilinear, uy Hed Lancclon Oc x 


Only; i.e. 
= = fe u, asc) (II-3) 


Tieegratangd Eouation.: (.i-3) by parts, 


R z= Ris RK e 
hj us dx = ie “if : x du,. (II-4) 
Assuming no slip at the wall, 
wetx/p c= 0: (II-5) 
2." 0 


From Equation, (Li-4) and (22-5), Equation (11-3) can be 


written, 
ay R 
aR = ce > du, (II 6) 
Or, 
eat R du, 
aR aa rz ok see) xs. (LT 7) 


The equation of motion (52b) for the laminar flow of Figure 
(II-2) takes the following simple form, 


VP = ="V% (ii=8) 


ae | 
= 


where P is pressure and t, stress tensor. Integrating 


Equation (II-8) for the present case, 
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eee fae ei BP zs 

( re AL x: (II-9) 
tT represents es From Equation (II-9) the wall shear 
stress, uae is obtained, 


AP 


Ti = ty R. (E310) 


Dividing), Equation (11=9) by (Ii=10) and solving for x, 


a fer) (II-11) 
W 


and, 


Sede (1112) 
Tw 


LE weesuipstituce (tik, +2) into (Li-7) , we get, 


qi 2 T du 
Peet eels Wa i hy 
3 
W 
OG, 
Qi R v Soe 
RR ge eee Glee deo eee 
W 


where u is average velocity and du, /dx is shear rate and 


will be designated by o, so that, 
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3 
se ae DIN Baksh gee (II-15) 
m = 2 0 
W 
For the assumed steady flow o is a function only of t. For 


Newtonian fluids, 
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bas ots rzserfa et xB\ ob bre yiinolsev opens a 7 we 
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z e 
(21-20 - 6 oth ee 


10% .7 30 Vigo notiszone? 


Zo 


Stee e (II-16) 
and for a fluid approximated by the power-law model, 
Ts nm, lo|", (II-17) 


where |o| is the absolute value of o, m, and n are rheological 


y 


parameters. By defining I as an average nominal Shear rate, 
Pee ee (II-18) 


Equation (II-15) can be written in the following form, 


t e 
tyes spel W to d_t. (T0119) 
a 0 


To modify Equation (II-19) for the multicylinder 


DOX , SL ESUSosucLicient , toguntroduce V/E 5 in place of Un! oe. 


Ss 


where E is surface porosity. Surface porosity is calculated 
as follows, 


F Ld = 
Se came aD Seeger 


Therefore, Equations (II-18) and (II-19) become, 


T =—=s-=- h S sou Ha (II-22) 


€ 


where the subscript e« indicates that the equation has been 
modified for the porous bed. Also, Equation (II-10) can be 
modified for the porous bed as follows, 


© ~P 2, 
Le Se K, (II-23) 
where L' is written instead of L, the length of the bed, 


to include tortuousity, l.e. 
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L' = Cy Ig’ GLI=24) 


Cy is geometric constant which accounts for tortuousity of 


the path of the fluid particles. The magnitude of Cy is 
estimated from Figure (II-l1) to be approximately equal to, 


Se yeas ya 
Cokie oees +0185) (II-25) 


where (1D/2) is half the circumference of the cylinder 

and N the number of rows. A more accurate value of Cy which 
will be presented later, should be determined experimentally. 
Substituting, (12-24) into (Li=23), 


PnP a 
ue = @L R. (II-26) 
4 
For a Newtonian fluid, Equation (II-22) by using Equation 


(II-16) becomes, 


ee : 
ee aul (II-27) 
solving for i, 
T 
oko (II-28) 


alt= 


The expression found for u is defined as Darcy's viscosity. 
pubstzttuting (Ii-22, 26) anto {LI-28), 


2 
Tesckaees © (ire 2o) 


ue ae Ty 
3Cy L 
By comparing Equation (II-29) with Darcy's law, Equation 
(I-4), the expression for the permeability of a Newtonian 


fluid in the parallel-plates capillary model bed is obtained 


as, 
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in which Cyr as stated earlier, has to be evaluated 


experimentally. 


II. 3 Power-Law Model 
Assuming fluid flow is purely viscous, the power- 
law model can be used to approximate the properties of the 


polymer solution. Solving for o from Equation (II-17), 


heh eee Gisteaa 
te 


where the minus sign is chosen because of the choice of 


coordinates. Substituting (II-31) into (II-22) and inte- 


grating, 
n Te Ln 
vee Teon arias ee 
ory 
1 + 22n n i 
be il a ea : (II-33) 
Further, we write iv and es in terms of K and C, using 
Equations. (li-—22, 26, 30), 
AP i 
TW SK/C\S ee (II-34) 
ee V/V3C 4 Ke. . (II-35) 


Tf we substitute (1134735) rnto (11-33) and solve” for Cae 


oe ‘eligder eine iseqeay AP (II-36) 
m @iaye ' 
dB ath 5 


Comparing Equation (II-36) with the modified Darcy's law, 


i.e. Equation (I-5), H can be obtained, 


20 


LU 


- yyy 
- ni! . ‘i 7 
ee We . . 
Ts Ma i 
He _ 
beasulsve od « ad 


ny a 
arone kv toga sot 


id 2© esistxweqorg eo otamleougge pp 


, {V5-TT) nettesps moak Ve) 2% ‘paivtos 
( Boe 
Yk 


tO. solods SsHs IO saueseosd nsagtS, ek opie aun 


“gjnr bas (Sie-EI) oval (1t-12) catsudttednl _ 
| aga * 5 


4 
i= 
~ - X 
sy ° 
oe eae ) . U 
- * — -s ~ bef 
4 wt be vad fa LO = tIa3 Ae = | 
4 - 
at a ‘ 
(cc } 


ne = ti 
.wel e'yousd beritibem sdt daiw (OREN) mi 


: as rteddo Pres nie a ee 


28 


1+2n 
n 


he afore sch amare (II-37) 
From (11-37), 2t ts evident that’ H =m = u-when-n- = 1. as 


for a Newtonian fluid. 


II. 4 Modified Friction Factor and 


Reynolds Number 
If we use Ergun's definition of the friction factor 


(54), we obtain, 
K adie * Seika (II-38) 


The Reynolds number is arbitrarily defined so that, 
f.. = 1/Npe- (ii 39) 


After straight forward algebraic manipulation, it follows 


Prommbquativons. Glie2, 30; 33) “and (1-5); 


2-n 
PEA aN. z 
Nre ROU We Hi : (II-40) 
4s 
For a Newtonian fluid, we set n = 1 and therefore, 
eee aes N's ‘ (11-41) 
Re 2 ae Cy E u 


II. 5 Application of Theory 

The permeability, K, of the bed can be determined 
from Equation (I-4) by a Newtonian fluid flow experiment. 
Then, from Equation (II-30) the tortuousity factor, Cy, can 
be predicted. Therefore, H from Equation (II-37), and con- 
sequently, AP from Equation (I-5) for the polymer solution 


runs are calculated. 


rus. peepat woitoty 
7 Q San 


tot Dz Ljokad oft Do ok thn ‘te @ “ous i ott » 


yi "i haiti 


; 4 


2s , t 7" - P ‘7 a 
‘ eo ae: Ob beALISeR Vili sigita te: @e su 


aewollot +i ,aolssalvqinemn sistiepls Beawror ee 
LV 
.(@-1) Bos. 4{8f ,0€ ‘SoTY) anoktal ~ 


. —_ 
yroanT 20 tiene tapes * rt a < 


ee i 


bornigiestebh sd nso bed sdd to .A Vat Lbdnomage a 


.jnemiteqxes wolt Stull apiaoswer «& vd (ser) not dst 
x tit 
MS ,,9 ,tostoet YILavoustos sis (Obit) snilaaa [ty ite 


a 


> 


; _ a 
“noo bas ,(T&-11) aotceupe moze i . saptesedt : «Be eee , 
we 9 


a 


Ws 
oliuioe xemylog ait soi (f=5) lad da most hoe 
| | tesatebien. 

® 7 a as 
7 


29 
Equation (I-5) predicts that a plot of AP v.s. V 

on log-log scale is a straight line with slope n in the 
absence of inertia and elastic effects. Upward deviation of 
experimental points from this straight line might be due to 
the inertia effect, elastic effect or both. Effect of 
elasticity would be confirmed by showing the absence of 
inertia effect with a Newtonian fluid of similar viscosity. 
An equivalent verification can be obtained by plotting fric- 
tion factor v.s. Reynolds number which according to (II-39) 


would give a straight line with slope -1l on log-log scale. 
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CHAPTER III 
EXPERIMENTAL EQUIPMENT AND PROCEDURE 


Pita eMulta-Cylinder Matrix 

To simulate flow in porous media a regularly spaced 
matrix of 6mm diameter glass rods in a triangular arrangement 
was used. Figure (III-1) illustrates a unit cell of trian- 
gular pitch. The cylinders were contained in a rectangular 
channel of inside cross section 5.2 cm x 6.5 cm--Figure 
(III-2). The top and side walls of the channel were made 
of 1/2" Plexiglass sheets and the bottom, inlet and outlet, 
of stainless steel. On the bottom plate a lmm brass sheet 
was set to keep the cylinders fixed in their positions. The 
inside surface of the rear side wall was painted black by 
"Wrought Iron Black" paint to eliminate reflection. The 
front side of the box (i.e. the side facing the light source) 
was covered with a black hard paper which had an approxi- 
mately Imm horizontal slit at the middle to provide a narrow 
slit of Might from the*broader slit ‘coming/from the light 
source. 

To minimize the inlet and outlet disturbances, two 
distributors, 26 cm apart, were used; one before the first 
row and the other after the last row. The distance between 
the distributors and the inlet/outlet plates was 1.5 cm. 

The length of the model was 40 rows of of cylinders, each 
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FIGURE lll-1 TRIANGULAR-PITCH ISOCELL 


FIGURE |II-2 MULTI-CYLINDER MATRIX 
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row alternatively 9 or 10 cylinders across. The minimum 
and maximum pore openings, shown as d and b in Figure III-l, 
were approximately 0.5 and lmm, respectively. The volumetric 
porosity, e€, of the bed was 0.357. Measurement of e€ was 
performed by measuring the void volumes, vo° The box was 


filled with water, then the volume of the Water, jie was 


measured. Knowing the total volume of the box, Vine porosity 
was obtained by, 
Vy Va 
Ee = Tm amie (00 Wa Os oa 
dk - 
Void volume can also be approximated by, 

Lt = Vin = Vo ’ (III=2) 
where Vo is the volume occupied by cylinders and is equal 
£0; 

"pine mD- 2 RAs, (III-3) 

¢ 4 


where Nin is the total number of cylinders. 

A similar model with cylinders of shorter length 
and larger diameter was used by Kyle, et al. (55). A smaller 
diameter was chosen for the present study to achieve larger 
elongational rates (i.e. du /9z in. Equation V-3) which is 
required to observe any significant elastic effect. 

The model was constructed to represent a uniform 
homogeneous bed, but it was found that some inhomogeneity 
in the microscopic scale existed. Specifically, the minimum 


opening, which was nominally to be 0.5mm, varied between 0.46 


to 0.787mm--excluding the openings next to the walls. These 
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openings are given in Table (C-1l1) for three rows. And 

the diameter of rods which were nominally 6mm circular 
cross-section varied from 5.9 to 6.2mm. The cylinders were 
not perfectly circular and the maximum to minimum diameter 
ratio was as large as 1.03. 

The twentieth, twenty first and twenty second rows-- 
numbered from the inlet--were chosen to represent the bed. 
The positions are numbered from 1 to 32 as shown, schemati- 
Gally, inyFigure (211-3). JA complete map of these rows, 
measured at a total magnification 100 was prepared to define 
the spacings of the cylinders. This was performed by pro- 
jecting the pictures taken of the different positions which 
covered the whole area of interest. 

An attempt was made to obtain velocity profiles 
at the minimum openings, where the maximum pressure drop 
and the maximum deformation rates occur. These positions 
are along the line A-A in Figure (III-3). However, due to 
the difficulties involved in photography--e.g. not having 
enough light at some center line positions--it was not 
possible to have photographs of usable quality at these 
positions for all the 32 openings. As a result of this 
the complete distribution of the flow across a row had to be 
determined indirectly. Thus some of the pictures had 
to be taken from the positions below or above the 
line A-A, say B-B. The center line openings 


(i.e. the ‘openings along the line A-A such as 
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CD in Figure III-3 are designated by d; and the openings 
above or below A-A are called NCL, implying "Not Center Line." 
The relative magnitude of NCL with respect to d is an indic- 
ation of the relative position of the line B-B with respect 
to the line A-A. In the case of NCL, the orthogonal component 
of the streaks length was taken; but for the center line 
openings this need not be considered, because all the streaks 
are perpendicular to A-A. For the positions 1, 1l, 22, and 
3 2 the "Not Center Line," NCL, was taken even farther than 
B-B, say C-C, and those streamlines, crossing C-C, which 
passed through these positions were chosen. 

Every alternate row of the matrix had two large 

openings next to walls. These openings were in the order 
of magnitude of the radius of cylinders, i.e. 3mm, while the 
minimum openings between two cylinders as was already men- 
tioned, were designed to be 0.5mm. These large openings 
caused by-pass channels near the walls. That is a larger 
fraction of fluid passed from the openings next to walls. 


The by-pass area is shown in Figure (III-3) by dashed-line. 


III. 2 Experimental Fluids 

A partially hydrolized Polyacrylamide product of 
Dow Chemical Co., trade name SEPARAN AP-273 was used as 
polymer. This white granular powder is soluble in water 
and has a molecular weight of approximately three million. 
Its solution in water shows effects of elasticity and the 


rheological properties of 0.2% (by weight) solution, which 
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was used for the entire experiment, are given in Appendix 
A. Domestic tap water was ee for preparation of polymer 
solutions and care was taken to avoid mechanical degredation 
of the polymer molecules by stirring the solution very gently 
during preparation. The solutions were protected from 
bacterial degredation by approximately lcc of Sodium Benzoate 
powder in 10°ce solution. The specific gravity of the polymer 
solutions, determined for every new solution, was 0.9965 gr/cc. 

To provide a wide range of friction-factor Reynolds 
number for Newtonian fluids, Glycerol solutions ranging from 
96 to 58% Glycerol by volume in tap water were used. 

Aluminum metal dust, manufactured by J. T. Baker 
Chemical Co., Lot No. 2347, was used to visualize the flow 
pattern. The volumetric concentration of Aluminum dust was 


approximately 2cc of the bulk dust in ee solution. 


iii. )s) Pumping 


Fluid flow at a desired rate was obtained by a 
ZENITH gear pump driven by a variable speed transmission. 
For low flow rates (i.e. 0.685 to 6.85 cc/sec) one pump was 
used but for higher flow rates, which were required to 
increase the range of flow rate for pressure drop measurement, 
two pumps were used. The pumps were calibrated by weighing 
the polymer solution collected over a known period of time. 
Since the Glycerol solutions (i.e. more than 85% Glycerol by 
volume) are highly hygroscopic a complete closed recycle 


system was used. For the polymer solution, to avoid 
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degradation of the polymer molecules, recycling was avoided. 
Rheological properties were measured before and after running 
and consistency was checked to assure absence of any 


degradation. 


III. 4 Streak Photography Technique 

Streak photography not only gives the streamlines 
but also the velocity profile can be established by measuring 
the length of the streaks of known duration. These were pro- 
vided by placing a chopping disk perpendicular to the light 
beam to interrupt the light for known duration of time. In 
the present work the streak photography using a photomicro- 
graphy technique was applied to obtain the local velocity 


profile. 


ltl. 4-1) Optical ssystem 


The main features of the optical system are shown in 

Figure (III-4) which are briefly reviewed below. 
a. Light Source 

The light source used in the present study was the 
same as the one used by Rollin (56) and Catania (57). It was 
a CHRISTIE Xenon arc-lamp which operated from a three phase 
CHRISTIE transformer of D.C. maximum output of 3300 Watts 
at 100 amps. The necessary optical alignment and adjust- 
ments of the lamphouse No. BSF 50 were made to achieve the 
maximum brightness and light uniformity. Immediately in 
front of the lamphouse lens an adjustable horizontal slit, 


then a converging lens were placed. Therefore a narrow beam 
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Fluid Outlet 


Chopper Disk 


Variable Speed 
Transmission And 


Light Beam ots 


Converging Lens 


Adjustable Slit 


Light Source 


FIGURE IIl-4 SCHEMATIC PLAN VIEW OF OPTICAL 
SYSTEM 
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39 
to a thickness of, approximately, 2mm, less than the slits on 
the chopper disk but broader than the slit on the box was 
produced. The light slit was set horizontal so that the 
plane of light beam was perpendicular to the rods. 

b. Chopper Disk 

Two transparent lucite disks were assembled on a 
variable speed transmission through a pulley connection 
with a ratio of diameter of disk pulley to transmission 
pulley 2.7:1. In order to interrupt the light and provide 
streaks of known duration, a paper wheel with 8 equal radial 
slits and spokes was placed between the lucite disks and 
fastened together. The chopper disk was placed perpendicular 
to the light beam with the axis of rotation in the plane of 
the light slit. 

The rotational speed of the motor was set in order 
to obtain reasonable length of streaks. In an accelerating 
flow long streaks are not representative of the local velocity; 
on the other hand, uncertainty in the length measurement of 
very short streaks is more than that of longer ones. It is 
evident that higher flow rates require higher rotational 
speed. The rotational speed of the chopper above 100 RPM 
was measured by a Stroboscope Type 1531-A STROBOTAC. Below 
100 RPM, speed was determined using a stop watch. 

After the experiment was finished a disagreement 
was noticed between the actual and experimental total flow 
rates. The actual total flow rate was measured by the 


formerly calibrated pump, while the experimental flow rate 
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was calculated from the area under the velocity profiles. 
This disagreement appeared only for the two high flow rates 
where the rotational speed of the chopper was above 100 RPM. 
By checking the stroboscope, it was found that the stroboscope 
was not calibrated and the measured RPM's were 10.4 per cent 
less than the actual RPM's which was indeed equal to the 
inconsistency of the flow rates. 
C. Photomicrography Assembly 

As illustrated in Figure (III-5), a Miranda 35mm 
camera was positioned at the top of a phototube of a Wild M5 
microscope. The phototube was perpendicular to the light and 
the camera film parallel. Magnification on the microscope 
was 12X, but the real magnification on the film was 4X. 
Thirty-five mm Tri-X film with ASA speed of 400 was used. 
The exposure times were determined after a set of trials and 
were approximately 15 to 20 seconds. The exposure times 
were different according to the brightness of the position. 
For the positions where the light was not sufficient the 
exposure time was extended even to one minute in order to 
increase the probability of having a bright particle passing 


that.posie.on. 


III. 4-2 Measurement of the Streak Length 

A Bell & Howell projector was used to project the 
negatives on a large sheet of graph paper. Total magnifi- 
cation of the streaks on the graph paper was 100 in all cases. 
The distance between the projector and graph paper was adjus- 


ted so as to obtain the final magnification of 100. The graph 
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42 
paper had 0.1" divisions and the length of the streaks 
were read directly from the graph paper at the total magni- 
fication of 100. Therefore all the reported lengths in 
Tables (E-1/°2)°are in units’ of 107° inches. 

The system was calibrated by projecting photographs 
of a precision steel ruler with 0.5mm division. 

To eliminate the error due to the half-shadow at 
the beginning and the end of every light period of the 
chopper's opening, all measurements of streak length were 
performed on head-to-head or tail-to-tail basis. That is 


fnePlengen™ a” woe@mecasumed=noc  ™)%, 


le ae sagt RE jana ee Te 


Some typical photographs are given for Glycerol and Separan 
in Figures (IV-5, 6). 

Two major sources of error which caused the data 
points to be scattered in Figures (IV-7 to 12) or (IV-13) are 
how accurately the streaks could be located with respect to 
the edge of the rods and the hazy beginning and end of the 
streaks. 

By observing several photographs, points on the 
circumference of the rods were located and the cross sections 
of the rods were drawn onto the graph paper. In pictures 
which did not provide sharp edges, uncertainty in determining 
the points on the circumference is in the range of a few 


tenths of an inch in the magnified scale. This uncertainty 
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43 
is, approximately, 10 per cent the nominal center line 
epening, «div=.0i5mm. Noting «that .the,slope.of ;the 
local velocity profiles near the walls is high, the uncer- 
tainty in the location of the streaks causes significant 
scattering. For instance, the velocity profile of position 
No. 13 of Separan in Figure (IV-13) is a good example. In 
this profile, two groups of data points which belonged to two 
different photographs are observed at the left side of the 
mid-point of the opening. By shifting the upper group 0.05mm 
to the right, the scattering could be minimized significantly. 
This effect is more severe near the cylinders where the 
Slope of the velocity profiles is high, but is zero at the 
mid-point of the opening where the slope is almost zero. 

The formerly mentioned second source of error is the 
measurement of the length of the streaks which is due to the 
hazy beginning and end of the streaks. This could reach 
two-tenths of an inch in the magnified scale for good pictures. 
The percentage contribution of this error is negligible for 
long streaks, say 3 inches in magnified scale, but is 
significant for the short ones which are near the cylinders; 
e.g. for a 0.4 inches long streak in magnified scale, the 
percentage of the error in measurement is 50 per cent. 
Position No. 13 of Glycerol in Figure IV-13 is a good 
example of this kind of scattering. For the pictures of 
poor quality, the hazy part of streaks was much larger than 


what was reported above. This could reach to even half an 
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44 
inch in magnified scale. Some of the data points of 
positions No. 15 and 16 of Separan in Figure (IV-13) belong 
to this kind of picture. 

The "Not Center Line" positions were chosen so that 
both of these two sources of error were minimized. That is, 
they were chosen along a line where sharp image of the 
cylinder allowed accurate positioning and also where sharp 


streaks images occurred. 


III. 5 Pressure Drop Measurement 

Mercury and also two colored manometer liquids, 
MERIAN No. D-2883 and D-8325 with specific gravity of 2.95 
and 1.75 respectively, were used to measure the pressure 
drop. Pressure taps were located on the rear side of the 
box, 22.2 cm apart, between row No. 3 and row No. 38. 
Pressure drop was measured for pure, 96% (by volume), 81% (by 
volume), 58% (by volume) Glycerol, tap water, and 0.2% (by 
weight) Separan solutions. The range of flow rates for the 


different fluids is given below, 


PIL, Range of Flow Rates 
Pure Glycerol 0.431 - 3.42 cc/sec 
96% Glycerol 0.835°- 17.06 cc/sec 
81% Glycerol 0.835 = 10576 cc/sec 
58% Glycerol 2 elo oP o4eoD. CC7 BEC 
tap water 10.76 = 34,305 CC/7 sec 
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CHAPTER IV 
RESULTS 


IV. 1 Pressure Drop 

Figure (IV-1) is a plot of the measured pressure 
drop v.s. superficial velocity for the several’ Newtonian 
fluids. In this plot the empirical correlation of Bergelin, 
et al. (58) for the flow across banks of cylinders has been 
shown. This empirical correlation is discussed in Appendix 
B. The data points of this plot, together with the calcul- 
ated permeability, are recorded in Table (B-1l1). The per- 
meability of the bed has been calculated for individual 
rans Using Dauycy’s law, i.e% Equation (I-42) .." The arith 
metic average permeability of the runs number 1 to 41 in 
Table (B-1) has been chosen to approximate the permeability 
of the bed, i.e. 


Roce asa al Das eine (IV-1) 


Figure (IV-2) shows the pressure drop measurement v.s. 
superficial velocity for polymer solution. In this plot the 
predicted curve of pressure drop by Equation (I-5) based on 
the parallel-plates capillary model is presented. The 
empirical correlation of Bergelin which gave a good prediction 
for Glycerol solutions was modified for the polymer solution 


and the predicted curve is also given in Figure (IV-2). 
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49 
Figure (IV-3) shows the behavior of friction factor v.s. 
Reynolds number for Glycerol and polymer solutions. The 
correlations used in this plot are Equations (II-38, 40, 
41). The Separan data points of this plot and also those 
from Figure (IV-2) are recorded in Table (B-2). The Glycerol 


data points of Figure (IV-3) are presented in Table (B-1). 


IV. 2 Streamlines 

Photograph No. 1 in Figure (IV-4) presents the 
shape of the streamlines of 0.2% polymer solution across a 
Single submerged cylinder. The diameter of the cylinder 
was 3mm and the exposed length, 5.2 cm. The photograph 
has been taken in a plane perpendicular to the cylinder at 
its mid-point. The approach or free stream velocity was 
0.437 cm/sec. Photographs No. 2, 3, and 4 in Figure (IV-5), 
taken from position No. 25, present typical streamlines of 
pure Glycerol through the bank of glass rods. The flow rates 
were 0.685, 1.717, and 3.42 cc/sec respectively. Photographs 
No. 5, 6, and 7 in Figure (IV-6) have been taken from the 
same position and the same flow rates as Glycerol but with 


the 0.2% polymer solution. 


IV. 3 Local Velocity Profiles 

Local velocity profiles have been studied for three 
rows, namely row numbers 20, 21 and 22, and were established 
for total flow rates of 0.685, 1.717, 3.42 and 6.84 cc/sec. 


It was noticed that the percentage of the total flow rate 
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58 
passing through each opening was independent of the total 
flow rate through the bed. In other words, flow distribution 
remained the same for all sets of data. This suggests that 
by choosing a proper factor, the four sets of local velocity 
profiles could be superimposed. The data points of different 
flow rates are shown in the same plot in the Figures (IV-7) 
to (IV-12) by different symbols. For the sake of brevity, 
this has been performed only for one position of each row, 
namely positions No. 2, 17, and 27, and the velocity profiles 


* 
of the other positions are not presented here. In these 


plots projected streak length as measured directly from 
graph paper at the total magnification of 100 are plotted 
v.s. the position between two cylinders. However, the scale 
given in the plots has been chosen so that the profiles 
present the dimensionless interstitial velocity, i.e. u/V, 
v.S. position. The procedure to obtain the scale factor is 
discussed in Appendix C. 

The total data points obtained in the present study 
was over 5000, but the complete set of data points are pre- 
sented only for row No. 21 and minimum flow rate, i.e. q = 
0.685 cc/sec and for one position of each row, namely 
positions 2; 17, sand 27. ‘Again, for the: sake of brevity, 
only some of the data points are presented for the other 
positions and flow rates. These data points are recorded 

“The complete set of plotted local velocity profiles 


is available from the Chemical Engineering Department, Univer- 
sity of Alberta, Edmonton 7, Alberta, Canada. 
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in Tables (E-l1, 2). 

The area under the local velocity profiles as well 
as the d's (i.e. center line openings) are given in Table 
(C-1). It is evident that the area under the velocity 
profile--in appropriate units--multiplied by the length 
of the cylinder gives the portion of the flow passing 
through that specific opening, no matter whether the 
velocity profile is chosen at d or NCL. This flow rate 
will be designated by qi 

It should be noted from Table (C-1) that some of 
the values of d's of Glycerol solution runs are different 
from those of the polymer solution runs, although they were 
from the same positions. The reason is that, after the runs 
with polymer solution were completed, the system was taken 
apart and washed. Therefore, those glass rods which were 
not perfectly tight in their positions did not return to 
the same place as before. Also, as was mentioned in Section 
(III-1), the glass rods were not perfectly uniform; therefore, 
their relative positions could cause an increase or decrease 
ind. This, of course, occurred only for a few positions. 

The maximum interstitial velocity occurs at the 
mid-point of d--assuming symmetry in local velocity profile 
which is not exactly so. Non-symmetry in experimental 
velocity profiles, which is evident from Figures (IV-7) to 
(IV-12), is due to nonhomogeneity of the bed. For the 
positions that the center line local velocity profile could 


not be obtained, the maximum velocities were calculated by 
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Equations (D-15,17). The calculated and also experi- 
mental values are presented ie Table (V-3). 

In Figure (IV-13) theoretical parabolic velocity 
profiles for the flow of Newtonian fluids through parallel 
plates are compared with experimental data points for three 
positions. To compare with experimental data points of 
polymer solution, the velocity profiles predicted from flow 
of a power-law fluid through parallel plates are also pre- 
sented in Figure (IV-13). The Equations of these profiles 


are derived in Appendix D. 
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FIGURE |V-13 CENTERLINE VELOCITY PROFILE 
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CHAPTER V 
DISCUSSION OF THE RESULTS 


V. 1 Pressure Drop 

The linearity of the pressure drop with superficial 
velocity is evident from Figure (IV-1) which suggests the 
applicability of Darcy's law, i.e. Equation (I-4). Except 
for the water data the empirical correlation of Bergelin, 
et al. (58) to predict the pressure drop across an ideal bank 
of tubes is in excellent agreement with the data. The calcul- 
ations and use of this equation are given in detail in Appendix 
B. The water data are significantly higher than the predicted 
line. Since this is simply a Newtonian fluid it is expected 
inertial effects cause the discrepancy. Although the value 
of Reynolds number at which inertial effects become important 
cannot be uniquely defined for different porous media (14b), 
the particular value for a given media can be determined by 
experiments over a sufficiently large range of flow rates. 
Thus the Newtonian data as plotted on Figure IV-3 show that 
inertial effects become significant at a Reynolds number 
(defined by Equation II-41) of order unity. 

Assuming the flow of the polymer solution is pure 
viscous, Equation (I-5) which is equivalent to (II-36) can be 
applied. That is the slope of AP v.s. V on log-log coordinates 
is equal to the power-law parameter, n. The pressure drop 
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measurement v.s. flow rate for the polymer solution in 
Figure IV-2 indicates the pseudoplastic nature of the polymer 
solution. However, the slope varied from approximately 0.35 
to 0.92 over the range of flow rates while the slope of the 
viscometric data ranges from 0.35 to 0.44 in Figure A-3. 
Note in Figure A-3 the best fit lines have been taken over 
two portions of the curve to define power-law parameters used 
in the calculations. This is of little consequence in that 
calculated parameters, such as Reynolds number are not very 
sensitive to the exact value of the slope. The large slope 
at high flow rates and the fact that the Reynolds number is 
Significantly below that at which inertial effects are expec- 
ted suggest elasticity effects are occurring. Thus predictions 
for the pressure drop of the polymer solution need to be 
limited to the lower flow rates. 

Different approaches were tried to predict the 
pressure drop of polymer solution in the bed. First the 
modified capillary model discussed in Chapter II, then, the 
empirical correlation of Bergelin, et al. (58) with apparent 
viscosity was applied. For estimation of apparent viscosity 
Equation (B-38) was assumed to define the shear rate--refer 
to Appendix B. To predict P from Equation (I-5), one needs 
to estimate the viscosity level parameter, H. H for the 
present model is given by Equation (II-37) in which Cy and K 
must be determined, experimentally, by a Newtonian fluid. [In 


Section (IV-1) K was reported to be equal to 1.3 x Ly a om 
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67 
and surface porosity for the present model is calculated 
from (II-21) to be, 

ES = 0.134. (v= 1) 
This value is the arithmetic average of two consecutive 
rows, one without by-pass channel and the other with two 
by-pass channels near the walls. The values of K and a 
were used in Equation (II-30) to calculate Cy: 


C,2#—_0.682. (V-2) 


4 

The calculated pressure drop together with H are 
recorded in Table (B-2) and the complete predicted curve of 
AP v.s. V is shown in Figure (IV-2). 

In this approach the openings of parallel plates 
are taken equal to the minimum center line openings, d. This 
gives the highest possible value for a characteristic shear 
rate, therefore, smallest value for viscosity. In view of 
Equation (I-5), which is equivalent to (I1I-36), this gives 
the lowest possible value for AP. In fact, in the bed, there 
is a range of values of viscosity which occur and, therefore, 
the bed would have a somewhat higher value of average vis- 
cosity, and consequently AP, than the minimum which comes 
from the moan This problem does not occur for Newton- 


ian; because, uw is constant. In terms of f. - N the low 


Re’ 
value of viscosity causes higher Nee which makes the data 
points of Separan systematically shift to the right of 
Newtonian in Figure (IV-3). 


To further show the possibility of elasticity effect, 
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Deborah number was calculated. Deborah number was defined 


in Section (I-3) and is rewritten in the following form, 


SUncioe re. tat (v-3) 


Referring to the geometry of the bed du, /dz can be estimated 
via maximum velocity at center line openings. According to 


Figure (V-1l), 


amr) 


Az is approximately 0.33 cm for the present model and Ue 
is the average value of the maximum velocities at center 
line openings for two consecutive rows. To be on the 


conservative side, the shear dependent relaxation time was 


used for O- which was read from Figure (A-4). This required 
shear rate to estimate O-- The value of shear rate occurring 
at the minimum opening is of the order of, 
Un 2u 
— — Lf oa 
Oe yore dL | (V-5) 
and, 
tn ee (v-6) 
ae 
Therefore, 
2c), 
o=—>. (V-7) 
es 


Recorded values of o in Table (V-1) are the arithmatic 
: ‘ 2 
average value of two consecutive rows. That 1s q,/d at 


a given flow in the bed has been averaged for opening in 
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the two rows. 

Calculated Deborah numbers are presented in Table 
(V-1). According to Marshall and Metzner (45), the magnitude 
of Deborah number is high enough to expect significant 
elasticity effect. As mentioned in Section (I-3), Marshall 
and Metzner revealed that appreciable influences of the fluid 
elasticity were observed when Deborah number, defined by 
Equation (I-13), reached 0.05-0.06 or, defined by Equation 
(V-3), reached 0.2-0.24. It is noted that Deborah number 
defined by Equation (V-3) is four times as large as that of 
Equation (I-13). 

In Figure (V-2), fNee was plotted v.s. Deborah 
number for polymer solution. The upward deviation of the 
data points at high flow rates, is similar to that observed 
by Marshall and Metzner (45), although the analysis is not 


sufficient to indicate when this might occur in the figure. 


V. 2 Streamlines 

The elasticity effect predicted by Ultman, et al. 
(42) in the flow past a submerged cylinder is obviously 
present. in. photograph No. 1 in Figure (IV-4). The early 
bulge of the upstream streamlines, the earlier return to 
undisturbed flow downstream of the rod and the non-symmetry 
of the flow are characteristic behavior of viscoelastic fluid 
flow past a submerged cylinder. The approach velocity which 


was calculated from the length of streaks was U=0.437 cm/sec. 
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FIGURE V-2 DEPENDENCE OF ELASTICITY EFFECT 
UPON THE DEBORAH NUMBER OF THE FLOW PROCESS 
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Reynolds number is defined as, 
CCF pr i] (V-9) 


where density, op = 1, diameter of the rod, D = 0.326 cm, and 
ee 42 poise is the limit of viscosity at low shear rate. 
Therefore, 


Re, = 20700339: (V4.0) 


also, Weissenberg number is calculated by Equation (I-11) 
in which to be consistent with Ultman, et al. (42), 0, = 
3.126 sec was estimated from Bueche's theory, i.e. Equation 


(A-6). Therefore, 


Ne = 4.19. 1s) 


To compare the situation prevailing in the submerged bank 
of rods with that of single submerged rod, the following 
justification is made. The average velocity at the center 


line opening is taken as the approach velocity, 


q 
tl 
GS 
ll 


ees (V-12) 
Reie 
s 
The flow rates studied varied from 0.685 to 6.85 cc/sec, 
therefore, approach velocities varied from U = 0.151 cm/sec 


to 1.51 cm/sec. The other parameters remained the same. 


Therefore, 
Re, = 0.00216 to 0.0216, (V-13) 
N Se Ooo? CO /eovn (V-14) 
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By comparing the Reynolds and Weissenberg numbers of bank 

of glass rods with those of single rod, it was felt that the 

range of flow rates studied was broad enough to cover any 

elasticity effects observable in the shape of the streamlines. 
The initial goal of photographic study was to seek 

if the bulge effect for the polymer solution would be inten- 

sified by the presence of the other rods. The behavior could 

be, hopefully, related to the increased microscopic sweep 

efficiency in polymer flooding. However, on observation 

of streamlines in the bed it was apparent that simple 

determination of the streamlines would not clearly show 

the major differences in flow field of the two fluids. 

Therefore, the direction of study was changed and focused 


on local velocity profiles and flow distribution. 


V. 3 Visual Observations 

Visual observations were limited to the positions 
along rows 19 to 23 for Glycerol and Separan. The flow 
rates studied “were q ="0.685;7 19717, 3.42 and 6.85 cc/sec. 
The observations at points where the flow was anomalous 
revealed some significant differences between Newtonian and 


polymer solution. The positions where anomalous behavior 
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occurred with the description of the peculiarity are 
given in Figures (V-3, 4). These positions are where no 
flow exists, i.e. dead zones, flow is slow compared to the 
normal neighborhood openings, or the flow is in the opposite 
direction of the normal flow direction. It is noticed from 
Figures (V-3, 4) that the number of anomalous points for 
Newtonian fluid are 15 while those from polymer solution 
are only 5. The observed difference, especially the number 
of dead zones or nearly dead zones, might be a good quali- 
tative explanation of the increased microscopic sweep 
efficiency in polymer flooding. These points were usually 
near the walls where by-pass flow existed, however, this 
Situation is not unexpectable in the oil reservoir fields. 
The above observed distinguished difference imply 
that the streamlines are not the same for Glycerol and 
Separan solutions. This is also revealed in section (V. 4) 


by observing the difference in flow distribution. 


V. 4 Local Velocity Profiles and 
Flow Distribution 


It is clearly observed from Figures (IV-7) to 
(IV-12) that the local velocity profiles were not the 


same in the different openings. That is the small 
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DESCRIPTION OF ANOMALOUS BEHAVIOR 


1. Very small flow in the opposite direction of the normal 
flow direction at q=0.685, 1.717, 3.42, and’ 6.85 ce7aec. 


Very small flow at q=0.685 cc/sec. 
Same as (1). 


2 
3 
4. Same as (1) 
5. Same as (1). 
6. 


Almost no flow at q=0.685, 1.717 cc/sec. and same as (1) 
at q=3.42, 6.84 cc/sec. 


7. Almost no flow at all four flow vates-. 
8. Same as (7). 

9. Same as (1). 

10. Very small flow at all four flow rates. 
ll. Same as (10). 

12. Same as (2). 

13. Same as (1). 

14. Same as (10). 

15. Same as (10). 
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DESCRIPTION OF ANOMALOUS BEHAVIOR 


Almost no flow at q=0.685 cc/sec. and very small flow at 
gG=l.717,..3.42), and: 6585: ec/ see; 


Very small flow in the opposite direction of the normal 
flow direction at q=0.685, 1.717,°3.42 and 6.85 cc/sec. 


Almost no flow at q=0.685, 1.717 cc/sec. and very small 
flow at gq=3.42, 6.85 cc/sec. 


Very small flow at q=0.685, 1.717, 3.42, and 6.85 cc/sec. 


Almost no flow at q=0.685, 1.717, and 3.42 cc/sec. and 
very small flow at q=6.85 cc/sec. 
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inhomogeneity of the bed caused significant differences 
in velocity profiles. Also, the local velocity profiles of 
glycerol show larger differences between positions than 
those of Separan. This shows that the nonhomogeneity 
has less effect on polymer solution than Newtonian. To 
emphasize the differences in flow distribution, the maximum 
velocities at the center line openings of rows 20, 21 and 22 
are plotted against position for Glycerol and Separan solu- 
tions in Figures (V-5) to (V-7). The maximum velocities 
of the positions which could not be determined directly--due 
to the difficulty involved in photography at the center line 
positions--were calculated by Equations (D-15, 17). The 
justification for this can be found in Table (V-2) where 
the predicted and experimental maximum velocities are 
compared. In all cases, except one or two, quantitative 
agreement was obtained. No usable pictures could be taken 
from position No. 3 of polymer solution runs--neither at 
center line nor at NCL. The flow rate erecuen this position 
was calculated by subtracting the sum of the flow rates 
passing the other openings of row No. 20 from the total flow 
rate, 

From Figures (V-5) to (V-7) it is clear that there is 
a greater variation in velocity at the center of the openings 
with glycerol than with Separan. It should be noticed that 


the value of Ue ley for position No. 3°4n; Figure: (V-5), as 
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TABLE V-2 


MAXIMUM VELOCITIES AT CENTER LINE OPENINGS 
WweeeAM. DIMENSIONLESS 


Position Glycerol Separan 

No. Cake. EXp. Gales Exp. 

1 IRE) i SL Q 

2 E565 U8). 5 14.2 14.66 

3 14.16 ? 5 42 ? 

2 12.3% ? LOST ? 

= TH. 44 ? 8.64 ? 

6 14.8 2 ae SAG 
i 10.56 LO. 28 8.14 8263 

8 2 Ob £2). 32 Ted $405 

9 16.9 mF; OFS 2 10...28 
10 13.4 13.24 129 ? 
int 16.93 ? Loney 2 
ie G.o2 6 57 dee 8.64 
is Lic Oe a6 S06 6.16 
14 LOsS6 2 La Sa ? 
15 LO.28 Tb rib eyes 2) SAAS 
16 2 suk A Aprre G57 Sac 
Ah 10,02 9.36 6.29 ore ts 
18 15.53 14.8 | 8.06 9.45 
7 1208 ? 6.82 ToaF 


20 6.9 i 6.02 7.4 
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TABLE V-2 (continued) 


Position Glycerol Separan 

No. Gale. EXp. Cale. EXp. 
PaaS 4.89 CUI EN S23 Bed 
22 24.5 ? 34.6 2 
23 Zu 2 25 2 
24 10 10 695 8.85 
25 Lonts> 16.24 8.85 Onis 7 
26 Ll236 LO cs 8 6.42 1@e2 
Pg | 10.4 LOs (Ase Paez 
28 ae Oe) LOe2s 6. 79 BL 
29 6.703 ? 8.24 co teh 
30 15.06 ? 8.64 Oakes 
Ser 1354 NEST Ad Ure, 2 
32 ie 5 f 17.4 . 
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FIGURE V-5~ MAXIMUM CENTER LINE VELOCITIES 
ROW NO.20 
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FIGURE V-6) MAXIMUM CENTER LINE VELOCITIES 
ROW NO.21 
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FIGURE V-7 > MAXIMUM CENTER LINE VELOCITIES 
ROW NO.22 
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mentioned before is uncertain. Also, as was mentioned 
before, the openings of some positions have changed from 
Glycerol runs to Separan which makes a detailed position 
by position comparison without accounting for the openings 
invalid. 

The three curves of maximum velocities were added 
together to find the final resultant. They are presented in 
Figures (V-8, 9). Again, same kind of behavior observed for 
the single row in Figures (V-5, 6, 7) is noticed. The 
resultant curve for Separan in Figure (V-9) is much smoother 
than that for Glycerol in Figure (V-8). To clarify the 
qualitative conclusion derivable from Figures (V-8, 9), let 
us assume a straight undisturbed line of fluid particles 
before row No. 20. The shape of this line of fluid particles 
after tee through rows No. 20, 21 and 22 would be, 
roughly, the resultant curve--assuming that the major dis- 
turbances were at the minimum center line openings. 

it as alsor shown from Figures (V—-5, cy 7) that 
the flow rate passing through positions 1, 12 and 22, which 
are next to the wall, are significantly larger for the case 
of Separan than those of Glycerol. The increase 
in by-pass flow should not only be attributed to the 


wall effect, but the relative magnitude of the openings next 


— 


heaolinem esw aes 
io'xt bepaeio eved @ 
sf ; 
noltieog Belisjieh ®& Benen ae 
ne h. 
ae ~~ 
tnege edt tok pr iiuenss v0 tte noe iteg 
} ; . aM _ a a i" 
; . a 7 = a. 
t ; ' i 
tl 7 7 
srew asitiooloyv ae cunt to coves wig 
_ | 7 oe —s 
) a tq ete youl bined Semare faait only atk? ea uwize 
~ ¢ me i = 
eviszdo roliveresd to Drtr Smee sbook as e 4 sith siadlinhs 
2 _beotson @& (T .38 ,¢-V¥) sokuglt! ak wees 
_ FQ 
me 
come foum eat (@-V) etgypit mt 8 ESTSE 10% rane 
vytixelo of .(B8~V) atuelt af fo caoyia x08 
~ (2 Y) astvplt mort eldaviued nkneitnae atl 
, 
i oF 
[2 j ult to snail Bedyyierone stpheise 6 on mee 
a 
* , ‘ 24 a .. 4 
iieq uit to eall eid? To sqade edt 08 Om shed . e102 
, 
ed Blow SS Bas LS (Of .cm awos Apwo ats paige 
tam edi tedy prinvee 1s-~-sV2Ho tases Lueox odd 
: 12 
ninege sil 4stnae> miaminim eds as ox 6 851 
| | 3 ~ 
r AN -V) sexvplt moxt awode’ ote 2 vat 
‘ - - 
w% Si .f aeanotsieog mrnrrs pateesq etst wo? 
ys 7 
‘ ‘ Ap a ey Pree, 
(} xo sxopuel yfdnsottinapte os acetal s ot txer 
lai 
esont wait .ioweoylo To ezods y a tage 
_ 
, : - 
eit ot basediststs ed vino tor Ste £2 t ee 16g 
- . % , 7 oat _ ~~ ~o 
ByK go aft Zo 13008 qd ,toetts 
- 7 : Py 
a oO 7 a eo 


Legend: 


pene Row No. 20 
—e— Row No. 21 


—w— Resultant 


| nsionl 
Bayt elie sionless 


hot eo wea iS Sloe ie leah) 1956 202 1 
Positions 


FIGURE V-8_ RESULTANT OF MAXIMUM CENTER LINE 
WEILOGITIES? OP THREET ROWS = GRYCEROL 


a) 
3 
—_ 

a) 


es 


ie 


2! a — eee ee 


i? 


— — a ae — i 
- ~ - —_ - 


Ms 
6 a _ 
o *< 
r 
i 
7 a 
e 
: 
i 
rs je ¢ > ; 
i | —— 
aan ' - rg - 
s e iy 
s | C eT B ary O31 ¢e # € ToT 


enoitized _ i. a 


i aa 
WML ASWAO MUMIAAM AG TWAT ARO | 
407 ida aed AMT Ac 


< 


Legend: 


----a--- Row No. 20 
—e— Row No.21 
ee ON GNe 


—w— Resultant 


/V,dimensionless 


Ma xX 


eer ee i 6b {aL G thy oe We Roun POI eet 
Positions 


FIGURE V-9_- RESULTANT OF MAXIMUM CENTER LINE 
VELOCITIES OF THREE ROWS - SEPARAN 


7 & 
af 
~ 
OS ov wor a 
nurs ~~ ee | : 
f<.0A woh —e— 
: 
S$. .o4 wok - == 
insiivai —— 
ee 
™~, 
™~ 
= a r . 
a : 
tE 6 1 8S a : ia 7 
a ee ee a ee 9 - a: - a 
| cecoree etree e ae a £E If 
encitizos aa aa 
- 
a 
1) SSID MUMIAAM FO TMAT.JUe3s, e- V 3 Al 


WAR AS A2 ~ 2WOA sda 10 23) 


iy 4 
_ 


90 


to the walls and also the effect of the other neighborhood 
openings. 

To eliminate the effect of by-pass flow, the 
following procedure was undertaken. Positions 1, 2, 10, 1l, 
27 oul ely 2c, 2o7 ole and 32 which were near ‘the walls: were 
ignored. The flow rates of the remaining openings of each 
row were added together. Then the percentage of this total 
flow rate passing through each opening of a row was calculated. 
These values will be referred to as "corrected fractional 
flow rates", f. This procedure makes it also possible to 
compare, quantitatively, the flow rates of the individual 
openings in different rows. Since every alternate row has 
two by-pass openings and more flow rate will pass from these 
openings, therefore, less flow would be available for the 
openings far from the walls comparing with the openings of 
those rows without by-pass channel. However, by applying 
the above procedure, all the remaining openings of the three 
rows are scaled up on the same basis. Corrected fractional 
flow rates are recorded in Table (V-3), but fractional flow 
rates without applying by-pass correction are presented in 
Table (C-l). 


An attempt was made to find the behavior of flow 
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93 
distribution with respect to pore openings. Equation (D-8) 
predicts that the flow rate varies with opening to the power 
(i 2n)/mar Sincegna, P,aiteis. ckpected that the effec of 
openings is aggravated for pseudoplastic fluid compared to 
Newtonian in which n= 1. In Figures (V-10) and (V-11) 
corrected fractional flow rates for Glycerol and Separan, 
respectively, are plotted v.s. dimensionless openings, i.e. 
d/d in! where oie 2 = 0.33mm is the minimum d from Table (C-1l). 
The line with slope 3 is the theoretical equation for a 
Newtonian fluid through parallel plates, i.e. Equation (D-1l1), 
and the line with slope 4.5 in Figure (v-ll) is presenting 
Equation (D-8) with n = 0.4. It is evident from Figures 
(V-10,11) that the effect of openings on polymer fractional 
flow rates is somewhat less than expected theoretically. 
Furthermore a line with slope even less than 3, probably, 
will fit the experimental data points of Separan. 

Since the openings are coupled, i.e. the effect of 
the neighborhood openings can not be eliminated; therefore, 
the data points of Figures®(V-10, 11) are expected to be 
scattered. However, fa is a measure of local permeability 
and portrays the combined effect of the neighborhood openings. 


Since the fractional flow rates of Glycerol and Separan are 


directly proportional to local permeability, therefore, 


it Hsyifelttthay a plot of £, and f. V.S. f. do/d in and 
fo do/ain! respectively, will be less scattered. These 


are presented in Figures (V-12, 13). The eye-fit line 


passing through the points of Figure (V-13) has less slope 
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than that of Figure (V-12) which confirms the former 


discussion. From Figure (V-12) one can conclude, 


O27 il 

fa % (f. da) (V-15) 
or, approximately, 

f£. va? (V-16) 

G 
And from Figure (V-13), 

a 

fo ‘V(f, dg) (V-17) 

o~,..substituting.for fo Eromn.(V=16),, 
2 
Ey Ss (V-18) 


That is the fractional flow rates of the Glycerol solution 
is proportional with the openings to the power three, while 
those of Separan solution is proportional with the openings 
to the power two. 

The fact that the effect of openings is suppressed 
due to the use of polymer solution may have major effects on 
the stability of a displacement front: one via conventional 
dispersion and the other via dispersion coefficient. 
According to Scheidegger (14c), the most conservative 
dependency of dispersion coefficient on interstitial 
velocity is, 


Grebe u V1.9) 


where Cy is dispersion coefficient, bi, a constant, and u, 
interstitial velocity. This implies that by adding some 
little amount of polymer to the water not only the nonuniform 


interstitial velocities will tend to be more uniform, but the 
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99 
effect of non-homogeneity of the bed on dispersion coef- 
ficient would, also, be damped out. 

The above experimental observation reveals that 
the flow distribution of polymer solution in the bed can 
not be predicted simply by power-law model and some other 
effects such as elasticity have a significant influence. 
The failure of the power-law model in prediction of the flow 
behavior of polymer solution is also shown by direct 
comparison of the experimental data with predicted velocity 
profiles. The experimental data and predicted velocity 
profiles by power-law were given for positions No. 13, 15 
and 16 in Figure (IV-13). As it is noticed from the Figure 
(IV-13), the experimental data of Glycerol solution is fairly 
close to the parabola obtained from flow of Newtonian fluid 
through parallel plates. In contrast, although not shown 
vividly by the particular profiles plotted, predicted power- 
law bropaies are significantly flatter than the polymer 
solution data. This feature is best shown by Table (V-2) 
where the calculated values of ee pi which are based on 
power-law fluid flow through parallel plates, are, without 


exception, less than the experimental values. 
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CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 


VI. 1 Conclusions 
The following are the main conclusions of this 
study. 

1. Nonhomogeneity of the bed affected the flow distri- 
bution significantly. The effect of polymer solution was to 
Suppress variations in velocity with pore opening to levels 
even below those for Newtonian fluids. 

2. The experimental velocity profiles for Glycerol 
solution were consistent with prediction for Newtonian fluid 
flow through parallel plates. On the other hand those of 
the polymer solution were much more rounded than the predic- 
ted velocity profiles for a power-law fluid. 

3. Consistent with earlier studies, increased pressure 
drop because of elasticity were observed. However, these 
occurred at flow rates in excess of those at which velocity 
Measurements were made. Therefore, the elasticity of the 
fluid affects the flow distribution before it affects the 
pressure drop. 

4. Qualitative visual observations revealed that the 
number of the dead zones or nearly dead zones and recircul- 


ation regions for polymer solution were reduced relative to 
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the Newtonian case. This implies the possibility of 


enhanced displacement efficiency in a microscopic sense. 


VI. 2 Recommendations 

1. Two-phase flow, i.e. displacing oil with Newtonian 
and also polymer solution in the present bed, could be use- 
ful in visualizing fingering which, hopefully, could be 
related to the present single phase behavior. It is recom- 
mended that the by-pass channels be eliminated in the future 
model design. 

2. A study of local velocity profiles and flow distri- 
bution should be undertaken using the same kind of model 


with more diversified pore size distributions. 


‘6 
vy : 
fOr” , Saat J 
: i vs ny a | 
. ‘ho viilidtasog = fc gmt etd C 


7 
, ” > teks 
.8erea older oxoim & i youeloiti< on 
; a = : iy a 7” 
ye ‘a 
: epee 35. nds b- wo BLASTS 
vy ae ~ 


ok ie we A. scic 


7 (oan 
) ‘we, °. 
‘ if : 
4 earns 
aeinoswen tjtw [lo on 


Bib 
Lind of Bake bad 3 sneeeg ods ab not suton xo amy to% 

ad Bluop Wintegad dod pat xepnt? | unit 
-mooex ai 31 .stoivared sesig pipate pre ree 


~ 


by 4 
siustut eft oi bedsnimile od @ feansito, ae 


Ishom to Batt omsa ade pate neketseban od 


-eaoisudiatath esie etog Rad 


105 


ca 


BIBLIOGRAPHY 


Sandiford, B. B. "Laboratory and Field Studies of 
Water Floods Using Polymer Solutions to Increase Oil 
Recoveries. 8... fet. Tech., (August; 1964), pe 917. 


"History of Petroleum Engineering." A Publication of 
the American Petroleum Inst., (1961), p. 850. 


Muskat, M. "Physical Principles of Oil Production." 
McGraw-HilieBook CovwInc.}; (1949)yupam648s 

Aronotsky, J.:S. “Mobility Ratio-Its Influence on 
Flood Patterns During Water Encroachment." Trans. 


ALME PROS (D952). p. “dS. 


Aronofsky, J. S., and Ramey, H. J., Jr. "Mobility Ratio-- 
Its Influence on Injection or production Histories in 
Five-Spot Water Flood." Trans., AIME, 207 (1956), 
pe20 5% 


Hurst, W. "Determination of Performance Curve in Five- 
Spot Water Flood." (Petroleum Engineer, 25; P.B.-40 
(Apr it. 195398, 


Craig pe@ick, dr. j,eGeffien, of. emMepcandmMorse, GKopAciwioil 
Recovery Performance of Pattern Gas or Water Injec- 
tion Operations from Model Tests." Trans., AIME, 
204, H1955p, ep. 7. 


Dyes, Save Be pecandle,.Bo.Hij4and Erickson, »R. A. “Oil 
Production After Breakthrough as Influenced by 
Mobil®tvye Ration! OTranseeAIME, O200, 81954) 50p.281 . 


Mahaffey, J. L., Rutherford, W. M., and Mathews, C. S. 
"Sweep Efficiency of Miscible Displacement in a 


Five-Spot." Soc. Pet. Eng. J., 6 (March 1966), 
pak (Oe 

Slobod, R. L., and Candle, B. HH.” "X-Ray Shadowgraph 
Studies of Areal Sweepout Efficiencies." Trans. 


AIME, tIbeti vo), Ds» 200. 


Slobod. Rucle » and, TpOmds aRewhas Geicect, OL Transverse 
Diffusion on Fingering in Miscible-Phase Displace- 
Mente SOC. Pet. Eng... (March 1963) 7. p. 9 


102 


oO uy 

1%, - 
i’ 
Lee 

; 

ee 


fto assent ot enoltule 


.vLe , {db0@L ,taupak) 


tc 


“4 


" noitoubord £L0 16: ‘selqtonts4 ‘Saesael ia". 
: (¢ 


A ae 
aonsultal eel-oisast whe tkaean ‘we ted (2 


wheelie yei-xX" .H .@ ,etbasd bas ed 


in ghvial & 
Li rT 10 
7 


‘ 
135198 LQ] sé oldie 
Q 


eatbuse piok ee 


ne Ow, sine ae 
ia ee mpeg 


4 \ teat rms toxte% 


> ei iP 


oe 


ef) ORL, 00 Laat 


" Jtaom iosotorst 19038W palzud aa Jt64 ~ bog 
el 8 eet) | gee a 


LAL 


i 


Pa 
im 
. 


idow" ..aG yeb: «Hy yooss paws ha ‘s. 


—_ * 
_ er e"@: : = 


H motsjouboug xo noLjpetal ao epiteu y +t 
TOs , SMT ans xT "hoo o£ a “meat 1 
| 808 : q 
iY anit 
3V 102 (Si0ri107,xAI Fo noLéecbexesee® ae : 
ES ,xsenipad awelortet “. heolt tetew sole 


4 (eeee fesaa) 


M .? <5 tiep Lalita ce veh ~ 
$64 4a aoremcotsell yxevooeh 

enex? “.acgeoT Jebom ii anoizts7 qO nok: 

| oT gia f@@Ol) tts 

tne t 

KH ,moedtoixzd pee 1-H .@ ,Sfbasd oe . 

ooeiurt 26 dpvoudtise 76 sagas ioldeoke si] 

(h2@L) , 10S .SMYA .enea? * olen vdbdie [OM 


wedjsM Bos ,..M .W' ae e a ee se 
jnemesaiqeld aldis aim Re vonetok Daa qos 
J 


[ dowsM) 0 ,.U .pnit «det .go8 Rabon 


| astoneintz2a tyoqgsew2 Isaxh to ge 
.2aS .q q 9 (S80L) sel TA 


footie” .A .a \eamod? sium? suk 
~eldioaimM al prttaspale no 2 | 
¢l dox6M -L «past sted. ee dm 


x 


ae , 


as. 


Db. 


14. 


55). 


biG. 


Fs 


eke 


1 Hie Ves 


20. 


Zi 


226 


PSS 


24. 


Aas 


103 


Gaudle, B. H., and White ~Ms. iS. «Production Potential 
Changes During Sweep-Out in a Five-Spot System." 
TRAN Stes, AIMEE 2h6.461959)>, p., 446). 


Collins, sR. Bae {Blow sof, Fluids Through Porous Materials." 
Reinhold Publishing Corporation, New York: (1961), 
p. 174. 


Scheidegger, A. E. "The Physics of Flow Through Porous 
Media." University of Toronto Press, (a) p. 69, 
(bP: O0nn(ey 4 259 €l963)e. 


Christopher, R. H., and Middleman, S. "Power-Law Flow 
Through a Packed Tube." I & EC Fundamentals, 4, 
No. 4 (November 1965), p. 422. 


Pirson, S. J. "Oil Reservoir Engineering." McGraw- 
Ha lab+Book«Cow INGay (a) spe. DOC yea )uep. 6560, u61958).. 


Welge, H. J. "A Simplified Method for Computing Oil 
Recovery by Gas or Water Drive." Trans., AIME, 
£95) 911952) pap ae. 


Jones, M. A. "Waterflood Mobility Control: A Case 
History. WJ. Pet. Tech. als .(Septenber 1966),,.p. L151. 


Mungan, N. “Improved Waterflooding Through Mobility 
Control." Paper prepared for the 19th Canadian 
Chemical Engineering Conference and 3rd Symposium 
on Catalysis, Edmonton, Alberta (October 1969). 


Mungan, N.; pSmithy BeiWs,uand Thompson; «J. 4.1L.) some 
Aspects of Polymer Floods." J. Pet. Tech., 18 
(September 1966), p. 1143. 


Stiles, W. E. "Use of Permeability Distribution in 
Water-Flood Calculations." Trans., AIME, 186 (1949), 


PPrAgie 


Dykstra, H., and Parsons, R. L. "The Prediction of 
Oil Recovery by Water Flood." Secondary Recovery of 
Oi.) ain ebhesUni.ted States, APL, New York: (L950))7) p7io0. 


Pye, D. J. "Improved Secondary Recovery by Control of 
Water Mob@litys "sid.ebet.wréech.» (August, L9I64)7.p. 911. 


Gogarty, W.B. "Mobility Control with Polymer Solutions." 
Sacs Pet. Eng. J. 97, 40une;196/), p- Lol. 


DaubenuwaD. ohiegeMenzte, uD Ese. Floweons Polymer Solutions 
Through Porous Media." J. Pet. Tech. (August 1967), 
Der L0Gps 


£0 


fabdnestod sok sha 


me : 

> > 1 

r in lie my; 
mm 


" mote Ye sogteewdt Hi 


aah 


7 i) [eis 36M v ee 
; auoyt 204 d 0% ie f 


Rilo TOF 


wolt wsad-xewot” 12 . temo. bbbM bese seat A 


2 .g .(8@OLl sevpwA) .dtoet .3ed .h "YALL EdOM 3 xotat 


-anotsilos tomy iog isiw LowtneD. ysitidon" am 


. (Lael) txex wth is 5 en rox 


aa : = “tp No 
1 devout? wolt to sokeys @ oft" s0p 
22 .q (s). (sees adaomet Bo, Cait bast ie pce 
| | . . ° (fel) as — aah: ay 0 ri i a. ; 


.b ,alsdeomsbarT oe * i *.odu't bexoed | 


“WEBIDOM ”". perlasemi pid nto VIBE a aa 
BCR AD OSS eh > ae ee DRE -0 


[iO pattugem> xo? borttemM bakit Lqoika ar at a 
,aMIA ,.ene87T “,ewi3xd tedsW, xo aed vd Yxevo 
Te «gd a (S8@n) 


22655 A :LoxsmoD yIlLidoM boofiaedaw” 
¢. todmsS 2) Lee «§3t eT tet ater wey 
yiilidoM dpvoxzit pathooltsaeseW Ber voxqmn* 
ett boxrsqe 2g x98qed " Loxtne 
msigogmye2 brut bas sonetetnodD patieenion= Lleokn off © i 
.(@9e@l usdostoo) edz od A sQd noma yeheyledad): . 
emoc" .1 .L ,monqmodT Bas ,.W oo hg hme Ex tap 
BL ,.doeT .36% .L “.ebool® ramelot to etosgas 
, EbLIL: .g , (eel ele 


» a a 4 p Pd ~* 
iaevebows CSC. BAS WOT 


mk soldgudiztalia tit idseme®t to eet!" .8 ” 42 sel 
Qf) O8L ,SMIA ,-anat? "-erolssinola® Roolt-se7s a 


‘es ry a i 7 of 


yao ie at 


to moitolhesd adit" ll I aoseua bas Hh ye ptatya 
YXeVOISH yrs SAoDee -hoolt rejewW yd. yxev of 
q «{0G@l) :Ato¥Y wet .IGA .ested2 bedtat edd rm. 


reciadick aaa ‘aa 
GOR devel ee) ya Vie. OonpsaH yusbnosss beyoxqmi" G'S, 8 yd 


da 
oe 


a ’ T 
¥ ss FE ater. 
“ALS So 


| 
2 


-40L .@ ,{VeeL ont) F vet -PAa 


{jufos stemylol to wolt" .@ a yebsagh yak 
VoOl sam oA) -ffpeT .3e% .t * -£1 50M avorog 


ipyore 
~ iy _ 


ZO’; 


CAB. 


Pee 


29. 


30% 


chi 


Sa. 


o3% 


34. 


Soe 


3.6% 


Dae 


104 


Ssverer, Gy ha, and@arough (Ali (0S. aM. Miprediction of 
Sweep Efficiency in Polymer Flooding." Prod. Monthly, 
32, No, 10 "(October 1968), p. 22. 


Lee, K. S., and Claridge, E. L. "Areal Sweep Efficiency 
of Pseudoplastic Fluids in a Five-Spot Hele-Shaw 
Model." .s0c., Per. Eng. J. (March 2£968):, ip. 52. 


Smith, F. W. "The Behavior of Partially Hydrolyzed 
Polyacrylamide Solutions in Porous Media." J. Pet. 
Tech. (February 19 ), p. 148. 


McCartney, J. A., and Sloat, B. "Polymer Reduce Risk 
in Waterflooding." Pet. Eng., (December 1970), 
pF =74% 


Lozanski, W. R., and Martin, I. "Taber South-Canada's 
First Polymer Flood." Paper prepared for the 21st 
Annual Technical Meeting of the Petroleum Society 
of CIM in Calgary, Alberta (May 1970). 


White; VJt°R" , *Géddard,VJI0VE .) Cand.- Phila@ips ,GHsamMs 
"Use of Polymers to Control Water Production in Oil 
Wells." J. Pet. Tech. (February 1973), p. 143. 


Jewett, R. L., and Schurz, G. F. "Polymer Flooding 
--A Current Appraisal." Paper prepared for the 44th 
Annual Fall Meeting of the Society of Petroleum 
Engineers of AIME, Denver, Colorado: September 28- 
October 1, 1969. 


Ustick) R85. andoual house ,“0.0D 9 Comparison’ of 
Polymer Flooding and Water-Flooding at Huntington 
Beach, Calitornia, J. Pet. Tech.,. 19) (september 
LOGO). pe LLOSS 


Jennings, Ri’ *RV;?Rogers,'d.H. ,.andiWest;, TT. U. 
"Pactors Influencing Mobility Control by Polymer 


Solutions.?+) 26Pet Tech! (March 1971)9% p.i aol: 
BULCIK/ E2= Ue "A Note on the Flow Behavior of Polyacryl- 
amide Solutions in Porous Media." Prod. Monthly, 29 


NOI 6s” (Otine 1965)59 pei 4. 


"Pseudo Dilatant Flow of Polyacrylamide 
Solutions in Porous Media." Prod. Monthly, 3l, 
No .© 89 (March: 4/967 )&,)4-p.ei27 . 


Burcik, E. J., and Ferrer, J. "The Mechanism of 
Pseudo Dilatant Flow." Prod. Monthly, 32, No. 3 
(March *196895)" pia. 


a 


tO nolvoibea “t* oun a 
$1oM i bore |” -priboost 


» Vid 


oH Jtoq@-evil wa = 
-q , (Baer aa 

bexyi nad, vileta= sat ball eee 
: j .BLSemM Revo bieu me anoiszuloe ei nd — 
| aie intel vf |e bi: 


Arh 


i2iA eoubed aemylod" wa teeta’ Pe _ 2 


fio mai moLtouboxrd +teaeW fo vannye of exeooy lod % 
+f BI oO 4 ‘ ETe. 2 aud: 2%) moet , ged oe he 
At Ba 


notenlsauH 3s pitbool't= 


,(OV8l xetimoed) ,.pmt dee! * onde 
i : ‘ong 


i" 


Pye: 


mnsS-djuoz yeder” VF ve koxeM Brus. re 4 lens: 
. edt tot Hersqerq weged rool’ Yomyto® sai! 3 
1 - eal an - : a aw eave 7 
6 meloxrted ed? Io pattoon, tsotadtoeT fauna 
-(OVGL wem) BI od LA Quselsd ¢ OM ote 
; | ihe. ¥ ee 
.M .E ,agillid® bas ..8 .% Btpbbod: Jang! 


paiboolt usmyio@d”’ .% .5D ,sindedg Base «0 
iy st Baexvseqei¢ +ze7sd ". feelergga sts wy Bm 
iusloxrtet to ytelooe siz 20 palzaek Lis® fem nA 
+ Yodmssqea :obstolodD ,1svMed ,aMTA Fo | 
“taek: yf ele So 


avodt ize Dre 4 .! 
atBW bas enlpeane * 
aimetqe’) @! ..doeat .te@ Vt BR knee ery: 
HOLE aQy! = : 
Cae 7 


tC .. ,teoW bas ,.8 .U ,Stepam 4.8 10 een cine 


«uM 4 


zwemyiod ya Lozgaod yt lide rma age Yar ax0sDs z° 

-LCE -q .(1V@L dowsM) .doe? .268 +b *.anottufog 

| ei a aed 

t Lod Fo ~oOLVi wees wolt arist mo escH A": A es , 
YicitneM .bora “.sibeM evotod ak enolisyle 2 obs ms 

| -bL .g ,(@0@L opt) © ou 

ebimeily+ SBYLOT to wort Insts Lia TEP ng 
(f6 .¥tdteoM .boxd =". sibeM suoxod at anoksuloe 

Der mas | VENOOE ‘doxe “ie 9 

7 +) _ ast ; 

~ 

20 ae pprencs~e edt” .% vaovear bas ..t 6 

E © vi y SC (Yt te tc 4 »boxr"d ; " Wolt prey re bes 

, f ink sia icha os 


y 


Seem 


Bo). 


40. 


41. 


42. 


Ase 


44, 


45. 


47. 


48. 


49. 


5.0; 


LOS 


Burcik, E. J. "What, Why and How of Polymers for 
Water floeding Jio)Pet.ciEng. #i(Auigustialio6s jp". 60 . 


, and Walrond, K. W. "Microgel in Polyacryla- 
mide Solutions and Its Role in Mobility Control." 
Prod. Monthly, 32, No. 9, (September 1968), p. 12. 


Lynch, E. J., and MacWilliams, D. C. "Mobility Control 
with Partially Hydrolyzed Polyacrylamide--A Reply 
toa Emil Burcdk. "b. J/.9 Pet.) Tech.’ 21) YOctober 1969), 
ar 24 Ain 


Patvon;s J ie., Codts).« Ki. HA.) and, Colegrave, :G. ‘T. 
"Prediction of Polymer Flood Performance." Soc. 
Pet. ‘eng. Un )wit(Marchre 07s) stp 2e726 


Ultman, J. S. "Viscoelastic Uniform Streaming Flow 
Past Submerged Obstacles." Ph.D. dissertation, 
University of Delaware, Newark: (1970). 


Harvey, A. H. "An Investigation of the Flow of Polymer 
Solutions Through Porous Media." Ph.D. dissertation, 
University of Oklahoma, Norman (1967). 


Sadowski, T. J., and Bird, R. B. "Non-Newtonian Flow 
Through’ Porous’ Media.” Trans. Soc.’ Rheol.,, 9° (2), 
(L9G) HeMpeA 24132 


Marshall, R. J., and Metzner, A. B. "Flow of Visco- 
elastic Fluids Through Porous Media." SPE paper 
1687, prepared for the Society of Petroleum Engineers 


Symposium on Mechanics of Rheologically Complex Fluids, 


Houston, Texas: December 15-16, 1966. 


Astarita, G. "Two Dimensionless Groups Relevant in the 
Analysis of Steady Flows of Viscoelastic Materials." 
I & EC Fundamentals, 6, No. 2 (May 1967), p. 257. 


Gaditondepi NY. Y. seandemMdiddileman, S. ¥ “Flow? oh Visco= 
elastic Fluids Through Porous Media." SPE Paper 
1685, prepared for the Society of Petroleum 
Engineers Symposium on Mechanics of Rheologically 
Complex Fluids, Houston, Texas: December 15-16, 
1966. 


Harrington, R. E., and Zimm, B. H. "Anomalous Plugging 
of Sintered Glass Filters by High Molecular Weight 
Polymers." J. Polymer Science, Part A-2, 6 (1968), 
pin294e 

Jones, W. M., and Maddock, J. L. "Flow of Viscoelastic 


Liquids: Comparison of Departures from Laminar Flow 
in Porous Beds and in Tubes." SPE Paper 1686, (cont.) 


tol «et . Sob ai Xo eniusioeM no muleogaye 


Jes) ,é8a!l sz98aes 348 pedu mt hae akoe 4 


a 
roe Oar Le 
* 308 enenrylot — ou Bry ye Os 
8 2a tRaer sBuR eae 


~siyarosys oF mi 


+ LOAS AON 3 or 
es my fuao fe 


foxstaod wilitdan” aol ‘a Pent one fae, 
vidat &--obimeivros¢ilt ‘4 4: paeroreyn. yilsioas 
(@oe1 sedodon) [s , doe? ef .b Matai’ ~ 


synupe too Rtn wh 4 (OI 809 Py ag ete . 
' pomanroase? Bool renvlot to: Hottoihen 
cy wa UW vet Ao38M) ‘yeh — im 24 


T Opciones 4s oor otteed eh OME ene ' ee 
Jyepe hb .G. dt ' aolosted® bepxuamdwe 2 bre 
| y { ,areumheg 30: ‘wiezey 


pe 


a oe 
} } f I IOV ] fi ‘ ‘ . H A . 
) r " f f 2 5, a aT ences. 


nlio Yo vitereviAag 
ad 

| i bite BAA ysl v tate ro 6a 
ee ".eibeM avowed Apu ; 
«EDS oq leet 


a aoe —2r Hs 5” kt » 2. iberdes 1a 
woxud dovesd? ehivolt obtesioa: 
igtoue alt wot beaagaxd , sol 


~ 


4i~di wacdmeoed: seme 1 MOFAUOH, | 
i 


f inoLenamrad ow?” oO ‘gobo BEA 
C c wold ybses3 to sivylan 0 
ye a 0 ,2ladremeiat aie 


| ight e 
: | ' a: oe . & ,_ OSin ee cee See K_yshaos: 


4oOt Anvosd? ob hal? ney: 


IsksoR ot 7 de hs rsqeiqg 4 
yi | Bfi2eM ao muleogmya eager: j 
ri-el 3 rT ,qotauoR ,ebiol® stn Lame of 


5 oat EW 


rt : Sy , ] mi? bas «..28..8 
ts 5 tr cll vel aust lift eesld bexesd 
a | ei) } t ne fh a al ,enrsrod "Ts y¥io® eb "2 


oizjsalsooetLV to welt” +k -L ,A9CBbaM Ene ‘i 
walt sranimed mow cama! taged te pacha s magi 


if 
7 is os 
; = 7 bie 


oli: 


vie 


Doi. 


54. 


She 


56: 


= 


5G. 


pas 


60. 


61. 


62. 


63. 


106 


prepared for the Society of Petroleum Engineers 
symposium on Mechanics of Rheologically Complex 
Fluids, Houston, Texas: December 15-16, 1966. 


Savins, J. G. "Non-Newtonian Flow Through Porous 
Medias’) ‘lay, BG, 6Ge NO. 10 (OCtOnpet a7 69)5 p. 1bo. 


Bird; ke 6.7 sstewarc,; Whe, and mightrtoot, EE. N: 
"Transport Phenomena." John Wiley & Sons, Inc. 
Lape Oem ai 94 yc) De. 197 (1966). 


McKelvey, J. M. "Polymer Processing." John Wiley & 
Sonsslwinegs», €1962))t,inp>t a 00; 


=Srgun, S. "Fluid Flow Through Packed Columns."Chem. 
BGs 6OG-, 40 ,.NO.n2, (952); De 89. 


Kyle, C. R., and Perrine, R. L. "An Experimental 
Model for Visual Studies of Turbulent Flow in Porous 
Mater detiGc uC ne O pC ile Diniurt oo (hODDUaLVanlog key | De B10 


Rollin, A. L. "Similarity Laws and Turbulence Struc- 
ture of Drag Reducing Fluids." Ph.D. dissertation, 
University of Alberta, Edmonton (1971). 


Catania, P. Ph.D. dissertation, under print, University 
of Alberta, Edmonton (1973). 


Bergelinm, ©. Pscabedds,.K: 0.5 and=heaghton;, MatDs 
"Heat Transfer and Fluid Friction During Flow Across 
Banks of Tubes." VII "By-passing Between Tube Bundle 
and Shell." Chem. Eng. Prog. Symposium Series, 55, 
Nol, 297 959) ep: 1 4p 


Ishihara, K. "Incompressible Viscous Flow Across Banks 
of Tubes at Low Reynolds Numbers." Ph.D. disserta- 
tion, Oklahoma State University (1971). 


Boger ,o\D.« Veg4 and» Raman Murthy, AsaV . "Normals Stress 
Measurement and Evaporation Effects on the Weissen- 
berg Rheogoniometer." Trans. Soc. Rheol., 13 (1969), 
Ds 205) 


The Weissenberg Rheogoniometer, Instruction Mannual, 
Model R. 18, Sagamo Controls limited, Section D, 


Theory. 
Seyer, F. A., and Metzner, A. B. "Turbulence Phenomena 
in Drag Reducing Systems." AICHE J., 15, No. 3 


(May 1969), p. 426. 


Middleman, S. "The Flow of High Polymers." Interscience 
Publishers (1968), p. 148. 


ad fi; ' 


- @upentpas mo loys} 


ra Lenio> eid L HIS pe 
geet {ai—cl 
growe? opvasnit 
é - = \a 4 : 8 * } . a 
i | . 
vi ‘a 4 dy sik . $ » ; . 
eo%3 a voll  aieG) a rtirge cy en SxT 
(hae) TOL ast Cah gee oq: CM wD Uae 
DF - Fe 
> Fi ne Le tes Sol , me 
' on bees ie seeing £09” a _ y¥9¥ 
ry Q 1Sd) 1evora? \aiet 5 eS ‘< a 
: y og . Fen? = won La xe 
{ \ 7 Yow) 
: ‘, ri 4. ‘ 
; ii” ae : Sted ean be ne or ke OLY 
; We oH >is Oe | @ Eve LV. 203 Lak =" “ 7 
es 
( ; os yr tise ® bet) me 7 mes: eng 
; ' ne 
: : ‘tit he ra BL J Aad arte «A 
cr 
£ Ji veal Vo : 
; 9i463 Yeoar si i ctg “Si veae rr a2 
: i ic. ela ds 1 0g vad lA to 
B&B ‘ Lod * Ms: oO aps 
| hes 2037 yi 
: a? + 2 a 
\ + = | an _ 
; tt ety OB ma2 
(ah .q y (ROL) @S, 
; s 
, s[diees ,emeonl” ye 
ars Bblonynl wed. 3e 
. | Seon S2608 GMOS 
1 ? hero emest Bin 
i ) i ‘ i Jt 20g Vu Bree 
} be: { % & <i » a? J ~ ik . . 185 xO E 
, Jen f atl , w:wtamoLnovosdat prsdaeaaie 
, , ; Lt eloruteo® ome seed .8f 7 — 
rated ted: it 2 A «se0Rdet Bris a 2 at ys¢ 
a .ai y.U SMOTA " emedave patoubes ry Gs : 


7 a7 Lot ify 8s 16 wool’ eit a 


ae | — nin 


a 
ly De 


APPENDIX A 
RHEOLOGICAL PROPERTIES OF THE SOLUTIONS 


A WEISSENBERG Rheogoniometer was used to measure 
viscosity of the Glycerol and polymer solutions. It was cali- 
brated with a standard oil. Also, normal forces were measured 
for polymer solution to estimate relaxation time of the solu- 
tions. To eliminate the surface evaporation of the solution 
at the edge of the platen a thin film of a light nonvolatile 
oil was used as recommended in reference (60). However, 
inconsistency in the different sets of the measurements for 
the same fluid forced us to change the procedure and diameter- 
correction method, which will be discussed later, was per- 
formed to correct for evaporation. 

Plots of shear stress and viscosity v.s. shear rate 
are given in Figure (A-2) and (A-3). Figure (A-2) includes 
standard oil, pure Glycerol and 96% Glycerol, while the visco- 
metric behavior of polymer solution is presented in Figure 
(A-3). The straight lines of shear stress v.s. shear rate with 
slope unity obtained for standard oil and Glycerol solution 
indicate the Newtonian behavior of these fluids. It is also 
checked by constant viscosity horizontal lines obtained for 
these fluids. The curve of 0.2% (by weight) polymer solution 
indicated pseudoplastic behavior in the range of shear rates 
studied. The flow behavior index varied from 0.8 to 0.346 for 


low shear rate region and from 0.346 to 0.443 for high shear 
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rate region, :The apparent viscosity varied from 0.2 to 42 
poise at the lowest shear rate. Table (A-2) represents the 
viscometric data and the results obtained from the Rheogonio- 


meter for all the fluids used in the present study. 


I. Rheogoniometer Parameters 


From (61), the shear stress is given by, 


3T 3.82T 
a amar a (A-1) 
ZA OSS dp 
where, 
tT = shear stress, dynes/cm* 
Aa ie a 
r = torque = S,At Ky 
Sy = torque transducer sensitivity, micron/volt 
At = movement of torsion head transducer, volts 
Ky = torsion bar constant in dyne cm per micron 
movement of the transducer, 
Roy = radius of the platen, cm. 
A541 = diameter of the platen cm. 
The shear rate is given by, 
_ 360 en 
where, 
-1 
o = shear rate, sec 
« = angle of cone, degrees 


t = 27/8, second/revolution 
B = angular rotation of the platen, rad./sec. 


The normal force is calculated by, 
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where, 
Pain Po> = primary normal stress difference, 
dyne/cm 
F = S5 OV, normal force, dyne, 
S5 = sensitivity of the normal force trans- 
ducer, gr/volt 
g = acceleration of gravity, om/sec* 
Vig T OUuLpULS OL transducer, volt. 
The normal force transducer was calibrated by placing 
analytic weights on the plate while rotating and the sensi- 
tivity was found to be S. = 0.2 gx/volt. The calibration 
curve is presented in Figure (A-1) and the data in Table 
(A=1) 
Apparent viscosity was calculated by, 
wes, (A-4) 


and relaxation time (62) by, 


sc Vilui@ cee (A-5) 


Sie Ze 0 
The plot of relaxation time v.s. shear rate is given in 
Figure (A-4). The maximum relaxation time was estimated by 


Bueche's theory (63), 
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Lh. = viscosity of the solvent, poise, 
M = solute molecular weight, g=-mole 
C = solute concentration, gr /om> 
R_ = gas constant, gr om? /sec* g-mole ok 


g 
le 


It 


absolute temperature, ok 
The above parameters for 0.2% (by weight) Separan AP-273 


solution are, 


Ve. eo 42 poise 

Ue es 0.013 poise 
M= 3 xX aye g-mole 
C = 0.2/100 = 0.002 gr/em? 

Ry = 8.35 x Mie gx cm? /sec” g-mole ok 
Te! 2.9 2 Oke 

Therefore, 
Oe oe 1 2G ec: (A-7) 


O 
For the entire set of runs two different platens 
but only one torsion bar was used. 
Runs with small platen, 
dpi = 5 cm 
te 2S oy = 3 


K, = 106.6 dynes cm/micron. 
Sy = 1.008 micron/volt. 
S5 #022) Or/vVoLe. 


Therefore, Equations (A-1l, 2, 3) becomes, 
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Runs with Large Platen, 


d ‘= 7.5 jcm 

Ot Se DO) ye ie 
Ky = 106.6 dyne cm/micron, 
Sy = 1.008 micron/volt. 


Therefore, Equations (A-1l, 2, 3) become, 


toi= 0.974 At, (A-11) 
Dee 

(oy ae r (A-12) 

Pal - Po» = 4,443 Vis (A-13) 


II. Evaporation Correction 

Care should be taken that the fluid covers the whole 
area of the platens. However, due to the evaporation at the 
edges of the platens, the area covered by the liquid would 
always be less than the platen's area. Therefore, the 
platen's diameter in former calculation should be replaced 
by the exact diameter of the portion of the platen which is 


covered with liquid. 
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TABLE A-1 


CALIBRATION OF NORMAL FORCE TRANSDUCER 


Run Mass Volt 
No. (gr. ) (at) XT) 
1 Jk 4.9 
2 | 15 24 

S 5 26 

4 10 49.4 

5 20 Oo 

6 30 150 

7 50 247 
TABLE A-2 


VISCOMETRIC MEASUREMENTS FOR STANDARD OIL, GLYCEROL, 
AND POLYMER SOLUTION 


1; Standard Oil 
VISC. No. S-600-68-lh was used as standard oil, which had 
ViSCOsaty Ou 19/4 <C.P. at 20eCy (and, L294 jCoP. at 25°C. 


The temperature of this experiment was 21°C. 
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TABLE A-2 (continued) 


Run Oy At 2 2 H 
No. sec Volt dyne/cm Cabs 
S 2.65 digs Wns) 5099 1922 
6 0.848 5 16.42 1934 
7 0.265 Li 6 5526 1982 
8 0.0848 07.55 1.808 2430 
2. BurecGhycerel, T = 2b°C 
1 8723 130 12637 1376 
2 18.4 260 253595 1377 
3 SRS) 80 78 1356 
4 1.464 20.8 20925 1382 
5 3.68 ordi 49.7 1350 
6 os i228 124.8 1350 
NOTE: Run No.'s 1-3 are for the pure Glycerol before 
being used, while runs 4-6 are for the same solution 
after being used in the system. 
3. 96% Glycerol, T= 21°C 
L 18.4 LS Lon 609 
2 4.62 29 Bene Gi. 
3 aa3s 45 43.8 596 
4 Lis6 i2a8 TOO 608 
a, of Glycerol, T= 20°C 
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TABLE A-2 (continued) 


Run oy At v 5 \ 
No. sec Volt dyne/cm CEPe 
3 L4c13 | L108 78.3 
4 Shan aie| 105 43.1 TET RERS 
5 89 165 Giee5 75.8 
6 141.3 260 106.6 158 
5.9 58e3Glyceroblbs T = 20°C 
ans 44.6 LOE 4.31 9.66 
Z 71 LY G9 G83 
3 we 26 £0.67 Deb 
4 Lies Al 67.3 9.47 
6. Domestic Tap Water 
1 bE Las Omn3 9 1.04 
Z LZ 2<9 1.148 1A AF 
3 phi) 4.5 1.844 e039 
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i) 224 eto 230 Leos 
6 359 922 Be PT B olf 8 ol 
7. , 0.2%,.Separan AP-273 Solution; Test No.1, T= 20°C 
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NOTE: The sample was taken from a freshly prepared solution 
before being used for the experiment and the period between 
the time of solution preparation and the time of performing 


viscometric measurements was 15 days. 
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Run Cite At Pe, u ¥t wb” 224 0 
No. sec Volt dyne/cm CEP Volt dyne/cm sec. 
9 232.4 Goe5 bee 32.4 169 815 0.0233 
10 29225 V5 82.4 28.2 204 984 Ue0ZOs 
LE 368 84.5 96.6 2032 240 1188 Gy OG 7 
E2 462 94 O75 2332 262 P28 0.01306 
53 575 105 120 2029 252 £599 0.0101 


NOTE: The sample was taken from the same solution as Test No. 
1 but after experimental runs and viscometric measurements 


were performed same day as Test No. l. 
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14 116 49 61.1 52.4. 100 525 0.037 
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Run coe At T 2 u Vt Bie, C) 
No. sec Volt dyne/cm Gx. Bs. Volt dyne/cm sec. 
16 184 59 73 36 40 100 52.5 0.01938 
bi 232). 4 67 8336 36 b53 803 0.0206 
18 292). 5 7-5 90:5 3:0). 935 ay 893 0.01686 
ae 368 82 LO 2..5 21.8 180 945 O.. 025.2 
Z0 462 90 Lae 24.3 220 L156 O;. Ou L3 
2k 575 100 124.8 Dy 262 1575 0.0958 
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5 Di 249 moder: Shi sd 576 - ~ - 

6 3.68 14.8 16.25 442 < 2 S 
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Run ome At 7 3 L Ve é iiatgee: ) 
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18 232).4 68 84.8 3.0409 aes 656 0.01664 
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a2 S75 LOY 26 PANE SS) 260 1365 0.00942 


Test No. 5; “T4="f20¢2C 


Run O 


No. sec + Vole dyne/c” CePy 
L 0.202925 0.49 T2220 4200 
2 0.0462 Lies 1.774 3840 
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NOTE: Tests No. 3, 4, and 5 were performed on the samples 
of the same solution which was different from that of tests 
No. 1 and 2. However, the sample of test No. 3 was one day 
old and was taken from the freshly prepared solution, while 
the samples of tests No. 4 and 5 were 10 and 17 days old, 
respectively, and were taken after the experimental runs. 
The results of the different tests indicate the lack of 
major degradation of the samples. 
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APPENDIX B 
PRESSURE DROP AND RELATED CORRELATIONS 


Bergelin, et al. (58) suggested the following. 
empirical correlation to predict pressure drop of the flow 


across an ideal bank of tubes, 


m (B-1) 


For isothermal flow the viscosity correction need not be 
considered; therefore, 


_ 4K, u WN 


AP B 
Ne a (B-2) 
216 S, D 
where, 
We = 9 q = mass flow rate, gr/sec. 
Sp = dx 2xN_ = minimum free flow area, om 
Ne = number of the minimum openings in one row, 
N = number of major restriction encountered in 


flow through the bank, 
B = indicates an ideal bundle of tubes, i.e. 


without by-pass flow. 


K = dimensionless constant which is determined by 
the following procedure, 

£ = ae 7 | (B-3) 

B Re, 


where f. is: fraction factor’ and Re, is tube bank Reynolds 
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number, defined as, 


DW, 
Re, = aes % (B-4) 
B 
In general, for tube banks (59), 
P,-l 
t= 3 26 8 
fur 7 pe meee Ep) re), 
a Py 


This is obtained (59) by integrating the Hagen-Poiseuille 
equation for a creeping flow between two tubes and using the 


tube bank friction factor defined as, 
Ye 
2 KPO Sp 


fB ein (B-6) 


In Equation (B-5), F(P, ) is equal to, 


PPE) = ‘ig dx : (B=7)) 


[P, - 1 eee 


PL is pitch ratio and is defined as, 
S 
a z 
where Si is transverse pitch as shown in Figure (III-1l1). It 
is worthwhile to mention that in staggered configuration the 
L Or Sin whichever is smaller) should 


be taken to calculate pitch ratio. For the present model, 


smaller pitch (i.e. S 


Sin ="O.Go CM: 3 (B-9) 
Therefore, 
Ea eto | (B-10) 


By graphical integration, F(P,) from Equation (B-7) was 


determined to be, 
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F(P,) = F(1.08) = 462.8. (B=Jab) 


Hence, by substituting into Equation (B-5), 
comparing (B-12) with (B-3), 
K, = L7i6n (B-13) 


As mentioned before, the Equation (B-1) was for an ideal 
bundle of tubes, i.e. without by-pass flow. To apply this 
equation to the present model, the by-pass flow should be 
eliminated. That is according to Figure (III-3), only the 
flow passing through the positions 13 to 20 should be con- 
Sidered. Hence, an imaginary wall is assumed to pass through 
the centers of the outermost cylinders which was shown by 
dashed-line in Figure (III-3). Therefore, ideal bundle 


volumetric flow rate, Gp is, 


dpi iSlBL04 poled + (B-14) 


i G =.0.556.q . (B-15) 


That is approximately 50% of the flow passes between the 
cylinders in row No. 21. For Separan, 


. 4-65 


dp = tay. Ob OLS 0.436 Qq e (B-16) 


In the above calculations, the flow rate passing through the 


imaginary wall and the outermost cylinders of row No. 21 was 
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neglected. A further assumption made is that the flow 
distribution would remain the same for the whole range of 
flow rates studied in pressure drop measurement. 

Considering the above assumptions, the Equation 
(B-2) can be used for the present model with the following 
known parameters, 

K, = 176 


If 
Wy Seni) O02 Gl 

N = 35 (i.e. the number of rows between the two 

pressure taps) 

S, = dxtxn,, = (0.05) (5.2) (9) = 2.34 om* 
D = 0.6 cm 
uw = 13.6 poise (pure Glycerol solution) 
u = 6.06 poise (96% Glycerol solution) 
uw = 0.765 poise (81% Glycerol solution) 


uw = 0.0964 poise (58% Glycerol solution) 


H 0.0104 poise (Domestic Tap Water) 
Substituting into Equation (B-2) and performing the necessary 


algebraic operations, 


AP, = 6.63 x toe q dyne/cm* (pure Glycerol) (B-17) 
AP = 2aeoo xX 10° q dyne/cm? (96% Glycerol) (B-18) 
AP, = 3.75 x 10°? q dyne/em* (81% Glycerol) (B-19) 
AP, Sage ex LO ag dyne/cm (58% Glycerol) (B-20) 
AP, = 50.0-¢ dyne/em* (Dom. Tap WwW.) *(B=21) 


In terms of superficial velocity, one can write, 


(B-22) 
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where A = 33.8 ene is the cross sectional area of the box, 


therefore, 
AP, = 2.24 x 10° V dyne/cm (pure Glycerol) (B-23) 
AP, = 1.0 x 10° v dyne/cm” (91% Glycerol) — (B-24) 
AP, = 1.27 x 10° V dyne/cm*® (81% Glycerol) —(B-25) 
AP, = 1.59 x 10° V dyne/cm” (58% Glycerol) —(B-26) 
AP,, = yo 8 la 10° V dyne/cm? (Dom. Tap W.) (B-27) 


Equations (B-23) to (B-27) together with the experimental 


data points were presented in Figure (IV-l). 


To apply Bergelin's correlation, i.e. Equation (B-2), 


to the polymer solution, an apparent viscosity must be defined 


for each flow rate, which the flow behavior indices, i.e. 


™ and n, can be considered constant. 


As a first approxi- 


mation, these parameters are estimated directly from the 


experimental curve of pressure drop 
in Figure (IV-2). According to the 
region of flow rates studied can be 


Region of flow rate 


cc/sec. 
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Bie oo) 6 ae eee 0. 
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v.S. superficial velocity 
figure (IV-2), the whole 


divided into four parts, 


n m5 
346 10 

443 ino 
632 Po Ss 
O17 Oe ey 


n is the slope of the eye-fit line through the data points in 


Figure (IV-2). These lines have not been shown. m, was cal- 


culated by assuming power-law 
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m4 = po ‘ (B-28) 


To estimate the apparent viscosity, the following procedure 
was undertaken. 

The wall shear rate for flow through parallel plates 
is given by, 


Pate LAG 
Gy Ban ta e)a72 Vane 


where oO, is wall shear rate, Vv, average velocity and d is 
shown in Figure (II-2). Equation (B-29) can be easily derived 
from equation of motion for power-law fluid flow through 
parallel plates (53). To modify Equation (B-29) for the 
present model, v should be replaced by average interstitial 


velocity, i.e. 
Vi Se Say (B-30) 


and d/2 should be replaced by volumetric hydraulic radius 
of the bed. Therefore, 


oon be V a 
a eR, (B-31) 


on the other hand, R. can be calculated as follows (52c), 


H 


volume of voids 


Ru - wetted surface 
volume of voids S ; (B-32) 
_ volume of bed > a 
wetted surface V 


volume of bed 
where s. represents the wetted surface per volume of bed. 
Specific surface, Ss! is defined as, 


_ wetted surface . (B-33) 
Ss ~ Solid volume 
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For the present bed, 
ee (TDL) Nw : 
(mD*2/4) N 


(B-34) 


Ol 


T 
in calculation of Sor the area of the surrounding walls has 


been neglected. One can also write, 


S, = S, (1-e). (B-35) 


From’ Equations: (B=-32, 34, 35), 


a. gD - (B-36) 
ae raga 


Substituting (B-22, 36) into (B-31), 


4(1l-e) 2n + 1 
=> ——_:_ g. (B=377) 
W Ae“D n 
Knowing the following parameters, 


te e388 cas 


O7e0u 


€ 
Dies 206 Cm: 
Equation (B-37) becomes, 


Z 2h. i 
oO. = 0.995 ane qe (B-38) 


The estimated oe from Equation (B-38) laid in the region of 
Dekel LL00s0 sec +. Over this range of shear rates Figure: 
A-4 shows n to be approximately .346 to .443. After the 
initial estimate of shear rate using (B-38) by n estimated 
from the pressure drop measurements, these new values of 
flow behavior indices were used in apparent viscosity 
calculation. Therefore, the calculated apparent viscosity, 
consequently, the predicted pressure drops were based on 


the viscometric measurement of the fluid. 
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The calculated values of wall shear rate, apparent 
viscosity, and pressure drop are recorded in Table (B-3) 
and the predicted curve of pressure drop v.s. superficial 


velocity by Equation (B-2) is presented in Figure (IV-2). 
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APPENDIX C 


CONSTRUCTION OF VELOCITY PROFILE FROM 
STREAK LENGTH 
As was mentioned in Section III-4-2, the length 
of the streaks and openings were taken from 0.l-inch division 
graph paper at total magnification of 100. Data points pre- 
sented in Appendix E and the local velocity profile plotted 
in Figure (IV-7 to 12) are based on this magnified scale. 
However, the scale given in the plots has been chosen so 
that the profiles present the dimensionless interstitial 
velocity, i.e. u/V v.s. position. The RPM's reported in 
Tables (E-l1, 2) are based on the chopper disk with one 
radial slit. That is, if the experimental RPM was 25 
revolution per minute with a 4-slit chopper disk, it would 
be reported 100 RPM. Because the time duration of every 
period of opening or closing the light would be same for a 
4-slit chopper at 25 RPM as for one-slit at 100 RPM. For 
every rotating disk, one can write, 


60 
RPM 


where t, is the duration of one revolution. For a one- 


t= sec/rev. (C-1) 


slit chopper disk, 
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the light. Therefore, taking s as the length of streak, 
the velocity is calculated by, 


tp 30 


(Cm) 


The flow rate passing through each opening is 


obtained by, 


erie =e Madey (C-4) 


where x-direction is shown in Figure (V-1). Substituting 


(C-3)} Tntco (C-4), 
= ——— 2f- sdx. (C-5) 


The area under the velocity profiles, ALs is obtained by, 


d a RPM f sdx.. (C-6) 
V fe) 3 O 


However, the plots in Figures (IV-7 to 12) are u/V v.s. 
x. Therefore, the area under the profiles, As, which have 
been presented for three positions in Figures (IV-7 to 12), 


is given by, 


_Ay_ RPM .d ‘ 
A. = V = 30V 0 sdx. (Ee 14) 
A.'S are recorded in Table (C-1). Comparing (C-5) with 


(C-7), q, can be obtained in terms of A,, 
She cee EN, (C=8) 


It is understood from Tables (E-l, 2) that all the lengths 
are reported in inches, i.e. 
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Kx = 25a x, 

Total volumetric flow rate can be obtained by using 
sum of the A,'S for each row in place of A; im Equation (C=-8). 
These values were in a good agreement with the flow rates 
measured by the formerly calibrated pump. For instance, 


for the Glycerol solution, 


row No. 20 7.> »} A. = 6.632 cm. 
row No. Zl S253 7), A. = 6.497 cm. 
FOWeINO. 22 Cope A. = 6.567 cm. 


Therefore, the calculated flow rate from Equation (C-8) in 
the case of minimum flow rate, i.e. q = 0.685 cc/sec. or 


V= (0.685/733.8) cm/sec. 


LOW NO. 20s? >> 6G 0.698 cc/sec. 


row No. 21 2-44 gq 0.684 cc/sec. 
row No. 22 >-2+ q = 0.691 cc/sec. 


and for the Separan solution, 


row No. 2b 34>> 2 A. 5 oe O° CRI. 


Gesb4aZ Cir. 


LOW NOwe co oo 1. As 
Therefore, the calculated flow rate, 
row No. 21 +4+>> q = 0.628 cc/sec. 


row No. 22 27->> g = 0.648 cc/sec. 
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TABLE C=-1 


OPENINGS, d, AREA UNDER THE DIMENSIONLESS VELOCITY 
PROFILES, Aj, AND FRACTIONAL FLOW RATES, q3% 


Position dio’, cm. A,x103 “i qs 
No. Glycerol Separan Glycerol Separan Glycerol Separan 
ROW NO. 20 
1 iD ie Lies 1891 18.4 314.83 
2 >.69 528. 720 562 List 9.29 
3 5.08 5.08 488 ? i Pesaylh 4 
4 5.74 408 493 400 eek) 6.61 
5 5.64 4.87 437 340 6.473 5638 
6 6.7 G.7 673 451 LO Sa 7.45 
7 5.49 5.49 392 387. 6.04 5557 
8 Serial DL 467 SZ 7 eo 5.16 
9 area 6.66 SOL 528 LES Wateie 8.74 
10 4.9 4.9 446 BZ 6.86 8.46 
A 3.68 3.68 422 440 6.2 YB 8 
ROW NO. 21 
1, 40 CW AATE L877 2380 2Bag S925 
23 ete) ee oRS) 3.2.0 249 503 4.9 
14 6.6 6.6 509 402 7.84 62902 
IAS, Ore Siegoy k 424 374 6.53 6. li 
16 65d j2 559 502 297 en 4.91 
7 S220 4.96 356 Psa Ps 5.49 ed Le 
18 6.76 6.76 veo 418 LOL 33 6.19 
19 S21 area 401 288 6.17 4.75 
20 5.44 5.44 Boz 286 oo 4.73 
7.6 2ea0 3240 1010 1010 5.56 16.64 
ROW NO. 22 
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TABLE C-1 (continued) 


— sss sss 


Position dx10°, cm. oe aie qs 
Nol Glycerol Separan Glycerol Separan Glycerol Separan 
23 oO; 92 5.14 991 1038 A ba ely dig open BG} 
24 5.59 bwo9 378 300 5.86 4.96 
25 637 6.7 790 461 2G LOL 
26 5.41 5.08 417 301 6.41 4.98 
20 otic oN) pe heaat 399 318 6.14 5.26 
28 SPaaas SL 369 298 5.109 4.26 
29 6.76 6.76 734 451 Lies 7.45 
30 5.46 5.46 Sou 381 Soy, O29 
ole 4.96 4.96 509 510 7.84 8.42 
ie! Sa 32'S 437 464 6.73 Thee 


NOTE: q% for Separan are based on average LA; of row No. 21 
and 22, while q% for Glycerol are based on the minimum total 


flow rate, i.e. q = 0.685 which is equivalent to A; = 6.49 cm. 
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APPENDIX D 
FLOW THROUGH PARALLEL PLATES-~-~POWER LAW MODEL 


The equation of motion can be, easily, solved for 
a power-law fluid flow through two parallel plates to obtain 
the velocity profile. The direction of the coordinates are 
given in Figure (II-2) 
rn 
tan —— 


i 
a n 
n 


M3 x 
==  -T (AP R n [1- (3) j 


Z 1 UA ORT) Oaks 
where m, and n are power-law model parameters and R = d/2 


shown in Figure (II-2). Also, 


a Le Silat © Oo) (D-2) 


After ingegrating with Equation (D-1) 


1 1+2n 
AB n (p=3) 


Equations (D-1) and (D-3) may approximate the rectilinear 
flow through the bundle of glass rods, assuming that the 
minimum center line opening between two contiguous rods of 
the same row, with length 2, is d and the flow rate passing 


through that opening is qa: Therefore, 
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where 


AG 
Equation (D-4) by (D-5). Hence, 
g 1l+n 
panda ne aL Cee 
Or en ee he aa 


where qe is known experimentally. 


profile can be established for each opening. 


setting x = 0, the maximum velocity 


can be obtained, 1.e. 
7 = o2n Wi 
max 1l+n Laas* 


From Equation (D-5) one can write, 


(1+2n) /n 


tie Crd 


where Ce is a constant. 
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is unknown, but can be eliminated by dividing 


(Dz6) 


Therefore, the velocity 


Also, by 


at center line opening 


(D7) 


(D-8) 


Equation (D-8) is based on the 


assumption that AP/AL is the same for all the openings. 


Equations (D-6, 7, 8) can be applied to a Newtonian 
fiuid’ by setting n*= 1, i.e. 
Q: 
ee eye 3 
Ue Ssvigat EL (Sailor (Deo) 
Se ie 
u => a 
max 2 la! (D-10) 
=C a (D-11) 
ne 6 : 
Equations (D-6, 7) can, also, be written in terms of the 
dimensionless velocity, 
fs q Ln 
Z. = bb2n *1 Oe 
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“max = it2n sia (D=13) 
V l+n Rdv- 
Using Equation (C-8), 
Lan 
u A. — 
+2 
A. 
+ 
mee = en a (D-15) 
For Newtonian fluids, 
u A. 
2 
y= 5 — 0-71, (D-16) 
“max 3 Ai 
= 5 4 (D-17) 
The calculated values of a A by Equations (D-15, 17) are 


recorded in Table (V-2) and the velocity profiles obtained 
from Equations (D-14, 16) have been shown in Figure (IV-13) 


for three positions. 


t-) vier : 
wv sdt bas (SV) 8 [dst ai 
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APPENDIX E 
EXPERIMENTAL DATA POINTS 


The complete data points of row No. 21 for the 
volumetric flow rate of 0.685 cc/sec, and the complete data 
points of one arbitrary position in each row will be pre- 
sented. For the other rows and other flow rates only a few 
necessary points to show the shape of the profiles will be 
given. The RPM reported in the following tables are based 
on the chopper disk with one radial slit--refer to Appendix 
C for more details. The pictures are labeled by the number 
of rolls and the number of exposure, E. The reported data 
in the following tables have already been converted to 
the same RPM. The necessary correction due to the error 
of the stroboscope reading mentioned in Section (III-4-1 b) 
is also included in the reported data. 

The length of the streaks, S, position X, and the 


openings, d were read directly from the 0.1 - inch division 


M! 
graph paper on which the pictures were projected and the 
total magnification was 100 in all cases. The procedure 
to analyze the data points was discussed in Appendix C. 

The position numbers are presented in Figure (III-3) and 


the direction of X and the origin with respect to the 


light and: flow direction isgiven in Figure (B-1) .. 
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Flow Direction 


Light Beam ) Cs ‘4 


FIGURE E-1 DIRECTION OF X AND THE POSITION 
OF ORIGIN WITH RESPECT TO THE LIGHT AND FLOW 
DIRECTION 
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xx10° 
Inch 


Sx10 
Inch 


Position No. 1 


ROLE No. 235 
NCL=91x1073in. 
q=0.685cc/sec. 


RPM=100 

El 
P2 Ged 
Pa nh! 
3h L239 
42 14.8 
47 ESS 
58 Lag 
61 Los 
82 14.8 
ou 14.8 

E2 
25 EE ok 
44 14.8 
51 Be Pee) 
64 14.8 
74 S45 
83 LSte5 


q=1.717cc/sec. 
RPM=250 


E6 
Ip 6.9 
34 LF d 
44 14.1 
54 15.2 
57 LoS 
66 LG? L 
76 151.95 
79 15'09 
of 14 

E7 
ne A 
30 ne ee. 
Dk LD 
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TABLE E-1l 


STREAK LENGTH DATA OF GLYCEROL 


Xx10 Sx10 
Inch Inch 
55 BS 25 
B3 54 
79 BS S1. 
86 14.6 
q=3.42cc/sec. 
RPM=500 
Ell 
10 6.4 
312 2:96 
3 3 35 
47 D4 22 
63 5 #4 
TA PS! e2 
86 2 
g=6.82cc/sec. 
RPM=1000 
E16 
77 LS) #6 
50 15.6 
Els? 
50 > 76 
55 15.6 
Sz deo te 
E18 
40 Miley "4 
5S 5 22 


Position No. -2 


Roll No. 239 
d=2.0).5x1073in. 
q=0.685cc/sec. 
RPM=100 


E16 
rete So: 83 
fhe Aa) 40.7 
E18 
12.4 44.4 


3 


Xx10 Sxi0~ Xx1L0 Sx10 
Inch Inch Inch Inch 
E19 NCL=62x10 ~in. 
7.5 AD.G BA58 5.8 
NCE=40x10-3in. 28 19.5 
E17 32 207 
10.5 18.5 44.6 t. 8 
22 DOW). Mae 5 
31 bie E12 
E18 3 sok 
25 2282 ACS 9 
E19 12.5 12.4 
16 22-20) TS Le 
E20 23.5 18.2 
32 1635 43.6 133 
NCL=62x1073in. 50 7.4 
E18 55 3.8 
26.7 19.6 E13 
pkg feat 20 12 He 4 
33.6 wees. 24 16 
35 Ish. 29 19.2 
E19 35.8 17 
19.5 48g . 3e°7 15.4 
Oe 18.5 44%3 12.6 
36.5 jna4 E14 
58.3 2.96 9 10-, 
25 16.7 
q=1.717cc/sec. BLS 
RPM=250 3.4 Pst 
d=20.5xl073in. 21.5 17.8 
Ell 50 9.6 
ial ih2. 26 
NCL=40xl0-2in! q=3.42¢éc/sec. 
Ell RPM=500 _4 
2.5 5.9 NCL=40x10 ~in. 
35.4 13 E6 
B12 7 Tieas 
eG Bar 38 a: 
6 10.1 NCL=62x1073in. 
E15 E6 
B87 G2 3.4 a 
v7, Tonk2  AepraA 8.5 
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TABLE E-l (continued) 


xx107 Sx10 
Inch Inch 
al TO7G 
44.4 10.9 
48 8.6 
52.2 5.4 
54.2 oer 
E18 
4.5 5.4 
15 15.4 
18.6 Te87 
a hd if) 
TEIG 19.8 
30.5 19.4 
39 16.8 
50.5 meg 
NCL=40x107-3in. 
Ell 
14.8 18.8 
24 24 
35 8.8 
E13 
4.6 10.9 
16 29°38 
El4 
ag 8.3 
34 Tos7 
E15 
tT? 8 16.8 
28 19.4 
34.4 1246 
E16 
AG 3.9 
8.5 15.2 
iby hee? 236 
31.3 ter? 
E17 
29.5 aD iy 
E18 
TPs D324 
14 yee 
083 pishe. 
E19 
a0.6 18.5 
E20 
16 yal Weg 
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Xx10 
Inch 


3 


RPM=1000 


Sx10 
Inch 
q=6.82cc/sec. 
NCL=62x10 °in. 
eee 
8.4 
E2 
el 
LQwS 
B hal ef 
Bie | 
E3 
DS 2 
Se: 
E4 
6.1 
6.7 
E5 
3 
Less 
So 
Sais: 
E6 
14.3 
5. 
6) 
E7 
4.4 
E8 
8.6 
16.9 
Ga 
E9 
‘ore 
15 
Ea s7 
F522 
oy hey, 
E10 
16.4 
6.2 
56 
3.4 
Now 259 
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Xx10 Sx10 
Inch Inch 
BL 
6 hey Ro we 3: 
a8 Ze 
E2 
9.9 9.9 

Meta’ 14 
42 BE i 
43.3 Bid eats: 
47 9.6 
E3 
pd Meg, 16.8 
Sa ec 
E5 
f) 6.9 
a Wada 8 Or: 
Ain p2 


Position No. 3 


Roll No. 263 
NCL=58x1073in. 
q=0.685cc/sec. 
RPM=100 


E16 
20 iy 
Be} 12.4 
33 Lag rai 
Sekees) sie 
Al 8 
47 ys 
so a peers 4 

EL7 
10.4 yk 
14.8 Ore) 
VE Weeds Died 
Chae: Ll 
OS) 1 fae 
36.8 la peg E 
a7. 4 6.6 
q=1.717cc/sec. 
RPM=250 

Ell 
3.4 aes 
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TABLE E-1 (continued) 


LE 


Xx10° 3 


Sx10 Sx10? xx10? sx10? xx107 Sx10 
Ech Inch Inch Inch Inch ncn Inch Inch 
8.8 5 ok q=0.685cc/sec. 15 9.9 
29.4 18.2 RPM=100 _., q=3.42cc/sec. 22 14.1 
39.4 9.4 NCL=57x10 ~in. RPM=500 38 Gis 
44.8 7 El E9 Ets 
48 5d 1532 Oy, 2.50 14.4 The 529 
53 oie 1 5 eat 3 11.7 44 Bret 9.8 dee 
E12 24.4 bMS E10 16 10.4 
5 Sv J S2e5 13,9 5 4-52 40.5 6.1 
rere} 4.7 Give Jet oe i 28 
bl 6.6 A3cc7 8 Ag .2 7.6 q=3.42cc/sec. 
27 hs 7 48 ~32°? 49.5 4.6 RPM=500 
38 O37 NCL=37x10 ~in. Ell E6 
47 579 10.8 15.4 9 5 17.4 Fal eo 
2975 4 c7 23 15.7 49 6.6 30.4 4229 
55 je, 30 Feed | E7 
B2 8.8 gq=6.82cc/sec. 5 4.2 
q=3.42cc/sec 34.5 5.5 RPM=100 26 14.3 
RPM=500 E2 E15 S45 Piss 
El 17.4 Love3 283 IS5 E8 
10.5 6.2 25.8 14.6 e264 LB 7 15 9.9 
Ee 10.3 28.3 11 24.6 P2S6 
35.5 Li 3 NGL=57x1073in. Positiongno. 5 
46.8 Se] 2525 1435 q=6.82cc/sec. 
E2 Sei. 2 La ROLL-Nox 265 RPM=1000 
8 4.2 44.5 8.8 gq=0.685cc/sec. El 
20 223 46 5.9 RPM=100 -3 29.4 14.8 
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Si. 5 8,8 E16 8.9 5.9 
q=1.717cc/sec 13.4 10 29 957 
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RPM=1000 mh JN 24.5 14 E4 
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ya | 2 LI.9. 4238 its 
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ey): 6a2 ll $.8 2065 3.3 
44 eZ 19 12.5 2% NG pe Roll No. 243 
E8 Cs yes 12.9 A3s5 4.4 q=0.685cc/sec. 
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E9 E6 q=1.717cc/sec. NCL=43x1L0 “in. 
5 B65 An! iS 7.2 RPM=250 E16 
3 79 Bie 6.6 Pies 
Position No. 4 24 1339 28 14,24 a ap a5 7 
34 13.6 42 Bcc SO] 22 
45 6.6 EZ 20 23.8 
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Inch Pnen Inca Inch Inch Inch Inch Inch 
oa) 2452 E14 5 252 
28 9 eo oA Ui piosge 8 2-26 Our eb 
32 ea EL5 L7 Leis Bo. ZOe3 
38 A i 14 CTs ae g bore 
E17 E16 2184 16.8 q=1.717cc/sec. 
8.6 226 oh eae Giz RPM=250 
34.8 E24 48 10.9 q=3.42cc/sec. mil 
RPM=500 1 6.8 
q=1.717cc/sec. q=6.82cc/sec. Bit geet Lo s3 
RPM=250 RPM=1000 Zee 340 6.1 2525 
Ell E2 5 83 244 292 
on E23 16 roel. 4G. Pers E12 
16). 6 24 E3 EG 35 LS e2 ee2 fe 
ES tai Ya pa 7 wero (2 ke 1532/5 3 NY ied 
Bu 16.3 A ing el 6 3.8 E12 3 ie 2325 
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20.6 23.4 EL E16 Ae 8 AGES 
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E7 q=1.717cc/sec. Position No. 8 E3 
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Bele 23 LOWE E16 7.3 3 eed 
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Bice 16.8 E7 9 10.8 Roll No. 246 
47.53 Wies 6 8.1 Bl? q=0.685cc/sec. 
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TABLE E-l (continued) 


xx107 Sx10° xx10° Sx10° xx102 sx10° xx10° Sx10° 
Inch Inch Inch Inch Inch Inch Inch Inch 
E20 90 Voie 2737 7235 294 27’. 
r,s 6.24 E3 E2 E7 
15 6.24 28 10.6 5.6 eZ 24 3.68 
43 13.22 139 pee? Age rr E8 
66.7 15.8 12.6 122: Tove 18.8 
80.6 TO.O2:. POSTLEIPONRTNO.| LS. lo .3 uid ea | E9 
aa omens | ANS. 13 52 4.4 ial gr 
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79 16.9 10.9 HOLT .2a.4 9.55 13.7 Led 
Er2 17-8 Leen 275 7.35 E13 
10 4.4 20.6 8.8 E4 9.55 Sigil 12.6 
16 5.7 E2 3 4.5 9 1638 
19 6.6 ve 15.4 9 ba NCL=32.9x1073in. 
31 Tror3 13.6 98) aie? 1228 Ell 
7 1655 2085 12.5 LOx 2 12.3 
q=3.42cc/sec. 23.8 2&9 236 7.6 i fsgeall PS 32 
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E6 Sine tase 29 5.9 6.5 To.5 
Te. 7 6.3 6.1 14.3 E5 Tra 1928 
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33 Albi gra A wg 10.3 8 nie ren 2722 6.6 
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29.5 10.5 21 Let d=25x107-3in. 
44 Toco E5 q=1.717cc/sec. E16 
68 15 2.6 6.25 RPM=250 4.8 9.6 
103 15.2 3.6 6.6 E6 NCL=70.8x1073in. 
124 ne +2 7 aed 6.1 7.3 8 1.6 
14.6 Lesa 32 13.6 t3.4 2 <8 
q=6.82cc/sec. 19.3 La s7 E7 20 4.8 
RPM=1000 NCL=32.9xl073in. 8 9.5 3o°2 7 
El El 29.5 221 B75 AD 
19.5 7.6 Bet 8.45 E8 64.4 22 
23 gee 8.8 es) as 13.6 E17 
32 10.9 12 foe Mo a. eee 24.3 5.5 
E2 14 15.1 d=25x1l0 ~in. 31 6.6 
20 gee 23.5 Pt E6 40.6 6.4 


,268\7956. = 
O02 
[ia 
i.vL C.€f 
fig 
a. St £.? 
8. aL e 
ak&-Oixe. Se=20K 
Lie 


ae? oat 


* 
.- « * 


£.0L 
£.8I 


bd 


*. » » . 


S 
é 
L 
LN 
8. 
b 
8 
é 
e, 
£ 
E 
t 
t 
an 
o 


é 
eg 
‘ 
“ee 
es 
r 
8 
it 
at 
es 


L535 


TABLE E-1l (continued) 


xx102 Sxl02 xx107 sx103 xx10° sx10? xx10° Sx10° 
Inch Inch Inch Inch Inch Inch Inch _Inch 
49 5.7 NCL=45.3x107in. 39 9.9 E4 
pA NE E13 19.6 16.5 
Position No. 14 6.8 : LE 35 27 Liss BS 
8.8 10.6. B80 B22 27 15.5 
ROLL No. 230 38 se? 3055 12.8 
q=0.685cc/sec. E12 q=6.82cc/sec. 
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10 2284 24.4 16.2 .39 AD VS <3 18 
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Das ites RPM=500 -3 35 5.6 5 6 
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13,8 22.6 22 10.9 d=24x10 “in 
E18 388.4 8.6 Position No.5 E8 
5 12.5 E8 225 1 
iy: 18.4 9.2 Sao (ROLL Nore 234 
18.5 24.6 28.5 15.2 qGe0.685cc/ sec. q=3.42cc/sec. 
E19 B05 12.6 RPM=100_., RPM=500 
iat POLL 41.6 9,3 d=24s10 “in: Roll. New Zod 
10 20.6 EQ E3 NELS42x10-3in} 
Lh. 6 26.9 5 S) 10.6 232 ET 
28 She Ps iS 36 $2) ..3 E4 2.4 3.2 
E20 2587 123.9 . DOs 26.5 10.4 PES 
5. % 16.2 By, 11.8 E5 26 13.8 
7.8 21 44 647 -. Le. 316.5 Boied 10.9 
45 225.0 47 5 NCL=42x10 “in. 34.4 7 
14.8 22a El 38.4 B. 7 
20.6 20.2 Roll No. 260 9 5 Pe E2 
23 LS NCL=59x1073in. 10.6 lhe A) 24.4 2.9 
24.3 if .8 Ell A ARE te 2 eg eit ee 4 Sak 
2265 O99: 620 14.3 39 Ze 
q=1.717cc/sec. 29 2 E2 
RPM=250_., 50 5.9 | 12 14.7 q=6 .82cc/sec. 
d=26x10 ~ ins E12 18 16.6 RPM=1000 
Font i 8 dS 1.18 8522 6.62 Ra 
25 Uez 19 8.9 E3 6 7 
4.8 15.5 20...6 10,6 12.5 14.7 26 16.8 
E14 27 286° F205 Aha th 33.4 6.4 
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TABLE E-l (continued) 


xx10? Sx10° xx10° Sx10° xx10? Sx10° xx102 Sx10° 
Inch Inch Inch Inch Inch Inch Inch Inch 
25 be ea 28 oe 25 14.7 
28 10.5 36.6 Ua Position No. 17 30 8.7 
E14 aes T285 E7 
30.4 G73 47 A.4> Roll No. 233 7 1s 
q=0.685cc/sec. 10.25 19,3 
Position No. 16 q=1.717cc/sec. RPM=100 3 20 14.7 
RPM=250 ag=19.5x10 “in. 27 Pt 
Roll Noys 23.2 Bld E2 E8 
q=0.685cc/sec. 20 £2.25 Zee Pied 1 eae es 2 
RPM=100 -3 2350 L553 5.5 20.6 20 14 
d=24.1x10 ~in 27 i248. i2.5 22.4 25 8.2 
Bay 30 42.9 413.5 22 E9 
20 ek 4 31.6 LB. 2 E4 8 
2 12. 38 Se as BS: 18.4 Wee a ee / 
E18 46.4 33-7 1gy 18.4 3 14 7 
25 6.6 a=24 1530 ~ in. E5 E1O 
E19 Bad 10 24.6 2.8 4.9 
2.5 oe a) 19.6 5 NCL=35.2xl0-21n. 4 6.6 
15.5 24% 2 E13 El 8 8.8 
ZL.3 OF oa 20:56 8.8 10.6 LOSS 2 2 
E20 E14 23 Ley: 20 12.56 
46.5 19.8 18.2 12.5 Zo 10.6 34 4.4 
N@L=52x107-3ant E15 E2 d=19.5x1073in¢ 
E16 ee) SA wae bl Ay .8 6.24 E6 
29 14.7 L9...8 LD / ooo Liss 10 8 24.2 
45.3 Ses: L.8 1 ie! E7 
Bi? q=3.42cc/sec. 22.5 14.7 he. 8 22.3 
eS 15.06 RPM=500 a 25'5 eL 11 20.5 
30.6 12 NECE=52x10 “ane§ 27.4 6.45 E9 
39 8.44 E7 E3 2 ddy 5 
45.3 5) 16.4 1.5 29 1538 
E18 22.15 14.5 E4 q=3.42cc/sec. 
oa .35 12.3 42 yD | bas 14 RPM=500 - 
58 12.8 48.5 2.3 26.5 11 NEbS=35.2x104 an. 
45 4.8 ow 5.5 EAL 
E19 q=6.82cc/sec. E5 14.6 12.9 
zZ0 13:6 RPM=1000 LOw3 2..s 16.4 1 Ae 
22 14.3 Be ES oe 15s k 20 +5.) 
26 15.4 25 Soe 25 12.5 
82.6 12 26 12 q=1.717cc/sec. Cw yo 4.2 
41 Seek 35 10.9 RPM=250 E12 
45.6 4.04 42 3/+8 E6 8 8.6 
E20 a=24 1x10 San 15 16.2 a3 i DE aa 
14 1302 Lb <5 Zz 19 15.4 16.4 14.8 
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TABLE E-1l (continued) 


Xx10° 


Inch 


Position No. 18 


Sx10° 
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Rolk No. 234 
q=0.685cc/sec. 
RPM=100 _, 
d=26.6x10 ~in. 
E16 
a5 20.4 
E18 
Zale E55 
8 29.6 
BS 302 
6.7 3226 
ZE.6 2057 
E19 
POS 3 eZ 
19 29 «5 
25.4 20.8 
E20 
5 2529 
3.5 32.6 
EOS 2 Sad 
Zu. 3 D2 
NCL=45x107-3in. 
E16 
Too i226 
18 23 
37 25 
E17 
3.8 EBS 
LS5 20.4 
ah 20.4 
29.8 20 
Ce, LS .2 
E18 
4.2 8 29 
Deed: 20 
Lees 219-8 
20 24... 1: 
Z2heo 20.4 
ie) Vers 
37 14.8 
E19 


xx107 sx102 
Inch Inch 
P22 E9.6 
18 pi: Doel 
2253 24.8 
23.8 24.8 
29 21.9 
3i2 18.9 
35.4 14.4 
40 925 
E20 
7.6 14.8 
9.5 16.3 
b2 Lous 
13145 Zod 
18.4 24.1 
20.8 2Ae6 
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q=1.717cc/sec. 
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Ell 
14.7 2202 
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xx10° sx10° 
Inch Inch 
RPM=100 F 
NCIS 44 .2x100° in. 
El 
27 20.3 
565 1285 
E2 
20.4 19.2 
25.6 13.5 
34 13.7 
40.3 5.9 
E3 
16.4 V6.6 
17.6 19.2 
20.7 16h6 
25 20.7 
34 14.8 


gq=1.717cc/sec. 
RPM=250 
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Shed Sie s 
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q=3.42cc/sec. 
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10 ive 
uO Dy Loss 
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Inch Inch 
2) 8.8 
q=6.82cc/sec. 
RPM=1000 
Roi” Now 227 
E2 
4 93 
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19.6 35D 
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it eee i 
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3G 
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E6 
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By 
43 18 
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1iL3 2.8 
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Xx102 sx10° 

Inch Inch 

NCh#@44 Sx102 in. 
E7 

28 es 
E9 

NG 12.2 


q=1.717cc/sec. 


RPM=250 
E4 
19 39-3 
NCL=58x10 ~in. 
E2 
oi Nbr) Sets) 
E3 
Mi? 5s 59 
Pee 8.6 
E4 
bag Gs: 7.4 
20.6 Tel 
S043 2...6 
45 4.9 
E5 
30 is ss 
CUS pes) SS 
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xx10° Sx107 

Inch Inch 
E12 

See D Te: 
BL5 

40 G.o 

Position No. 6 

Roll No. 210 

q=0.685cc/sec. 

RPM=100 _, 

NCL=43x1I0 “Gin. 
BL 

9 9.5 

¥.8 Nets: P26 
E13 

10 NORE 

25S Ya s7 
E14 

Lass Zc 

Zl, B47 
cl ie 

Bes: 9.2 

q=1.717cc/sec. 

RPM=250 
E16 

ZZ 14.5 

340 LZ atO 


q=3.42cc/sec. 


RPM=500 
El 

6.6 7.4 
E2 

17 1203 

26 15.5 
E4 

21aG 16.2 
E5 

7 een 

d=26.4x10 “in. 
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Ia? Pe hg?) 


xx10° 


Inch 


Sx10 


q=6.82cc/sec. 


RPM=1000 _, 
NCL=43x10 “in. 
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322 10.6 
E8 
Sed Be 
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near EO. 2 
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Sls iid ay 
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E7 
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4.6 L629 
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16 2ORe 
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6 16..9 
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oa 15 
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RPM=5 00 
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Pas) Lon Z 
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Position No. 8 


ROLL NOs c2dke 
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xx107 Sx10° xx10? Sx10° xx107 Sx10? xx10? sx10° 
Inch Inch tneh inch ‘ineh Eneh Inch Ineh 
Da 1S..6 B15 E19 
14 ES ..4 ley, 21.9 Positi6én No. 10 r 8.4 
E12 27, 16 
2 5.6 d=26.2x10-2inj4 R611 No. 214 q=1.717cc/sec. 
BLL q=0.685cc/sec. RPM=250 
q=6.82cc/sec. | 16.3 RPM=100 3 EY. 
RPM=1000 7 2026 NG@G=51L.52x10) cin. 0 8.9 
Ei BZ El 14 oa 
5 B79 234 983) 31765 139 47 5 
Liss 24 E2 HS 
Position No. 9 LOE) 32 3.9 0 LORS 
i3 23 PSU'S Arak E14 
ROLD No .':253 B53 20 Sree 
q=0.685cc/sec. q=3.42cc/sec. 36 le BLS 
RPM=100 ~3 RPM=500 39 Aes 45 4.8 
€d226.2x10 Zan? E6 45%2 6.9 
E17 6 22.9 E4 q=3.42cc/sec. 
5 20 ac E8 10 Ts RPM=500 
14.4 Zhe aa 10.6 48 $8.29 E6 
E18 L322 24.3 0 9.4 
16.2 27 2a NCL=34x1073in. Gel. 71/cc/see. 7 9 
E20 E6 RPM=250 Zo VT «6 
£4.23 2623 8 PS 4 E6 E7 
NCL=34x1073in. 26 HO.2 28 TL Gee DS Sa | 
E17 E8 E7 E8 
4.8 lee b 2 lien - 20 Li 47 4.8 
19 Eo.3 DTZ ONO) 22a a ees E9 
24.8 Boe 3 46 6.6 52 3.9 
E18 q=6.82cc/sec. E1O E1O 
19.4 OX RM. RPM=1000 23 14.9 43 4.8 
Zz 1832 Lgl 
E19 G25 V6.2 Position Nose q=6.82cc/sec. 
10.9 19 Sree L732 RPM=1000 
E3 Roll No. 215 El 
641.717 cc/sect aA 15.2 a=0.685 0 9.3 
RPM=250 Lae 16 RPM=100 Ee 10 B25 
Ell PAS) 82 10.2 NGL=68x10)~ in. E2 
5 bees E5 E16 20 6.6 
S.0 L5e¢ 135.2 i ee Ar a E4 
E12 d=26.2x1073in. E17 24 6.8 
3 8.8 E3 Hf) Ge. E5 
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TABLE E-2 (continued) 
xx10? sx10° xx102 sx10° xx10° sx10° xx10° Sx10 
Inch Tnch inch Tneh Aneh Inch Inch Inch 
95 20 102 19.8 E6 
Roll No. 138 E14 E2 64 of ee 
q=0.685cc/sec. 2 5eao 14.339 949 LS ak 140 i Peps 
RPM=100 _, 25 Te 2k4 445 20.3 E7 
d=146x10 “in. 3955 17.6 74 22 104 Ao Pes 
El 56 cys eae £3 E8 
15<5 12.95 "¥o4 L652 13 1t.9 38 19.8 
oi GS L513 E15 fod, 2 16.9 60 22.5 
32 i by Beet sta 25d BP a ie Sires) E4 90.5 23 
34 Lie Do 36 Lon 2s 18 PhS 96.5 20.9 
43.5 19.04 69 E96 Roll Now 177 
49 19.04 E16 Gar q=3.42cc/sec. 
59.5 2.0 tl S000 EO. 39 19 RPM=500 
82 20.95 45 17.3 60 21 E14 
oS 20.59 88 19.6 7H 21 48 214 
iMac we Lona? R07 15.8 108 16.7 ol 23.) 
19 15 SIGS PAS ie: 1d 130 LORS EL5 
E2 42 ead E2 2 24.4 
ee 13.04 WDD 20 43 Oe 2 E16 
1G 14.1 114 Dass.S 74 22 Ra 19 
28 a Je213 E18 98 19.25 30 20. 8 
55 20.4 35 16.95 127 55 62 Zens 
67 20 503 63 19.3 E3 EL7 
94 200 dy 81 2004 28 L714 32 20.5 
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E4 100 BS 72 Zee a 98 2ece 
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40 Sta 28 L6.2 53 Zk Rol No. 137 
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Aces 4962) “19 9.5 116.5 15 AG 20.5 
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64 20 El E6 Sh A030 
84 20.95 7 pa 209 ce LS 40 20.5 
8925 20 88 22 22 16.4 E20 
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%x107 Sx10° xx10°> Sx10° xx10° Sx10° xxl0” Sx10 
Inch Inch Inch mich ‘Inch Inch BE gvonal Inch 
43 ed 2 Roll No. 160 18 LOX el hora ies 
Ba ES PAN 16 
Positron No. 13°87. ig a RR In BOT poten E7 
Z0e6 a 23 Ee Ao 7.6 5 aS 
Roll No. 140 E2 Pil 2.50 4 8 
q=0.685cc/sec. Aa Oro E4 4D) 4 
RPM=100_, 10 ae wes ona E8 
d=25x10 ~in. 12°.°6 Wat De 5 1 fo 2 5.4 
E14 ED Pye nS 14.8 
ZU 6 A el eae" q=3.42cc/sec. Oe i lap? arg Zt EO 
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El Roll No. 159 ibpargs Vio Nd=6es2ce/ sec. 
ott ae] El ees Low RPM=1000 
5.4 14.7 di ehal 9 yep) E7 fae 
E2 Ste Ge 248 PA 18 PRs pee 
5.6 shi Ree E2 23.4 1A BRAY ey 2 
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E16 RPM=100 -3 a bay’ 263 5 tO 
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Xx10 Sx10 
Inch Inch 
BL 
9.5 os ae 
E3 
18 14.7 
6 TOME 
E7 
2445 ab 23 
E8 
3.9 8.4 
EY 
2.5 Be ae 
q=1.717cc/sec. 

RPM=250 
pcp aD 
a1 B22 
Bb2 
23 E35 
E17 
a. 5 ee 
19 16 
E20 
Lk La 
d=24x10 ~in. 
El2 
16 19 
E18 
18 16.9 
q=3.42cc/sec. 
RPM=500 
Roll No. 163 
Pole 
8 Grew 
E16 
ee ie LS. 7 
NCL=35x1073in. 
Ell 
Age D pe 
xe 4 
i. 5D 2.2 
F115 
GieZ ARS ew 


xXx10? 3 


Sx10 
Inch Dnach 
10 An 3S 
q=6.82cc/sec. 
RPM=1000 
E7 
S25 9.8 
E8 
4.8 cheat 
8 i Nas 
E10 
23 a4 13.4 
d=24x10 ~in 
E8 
eee 1G 
Position No. 16 
Roll, Nos 7146 
q=0.685cc/sec. 
RPM=100_, 
d=22x10 “in. 
El 
528 L226 
E2 
13 Tek 6 
E3 
6.3 ibe 
13 18.6 
E8 
12 19.6 
1S)..5 UT 
E10 
Os 8.16 
Roll No. 166 
BLL 
PyAw i eee 
E20 
3 222 
NCL=52x107-3in. 
Reb 
4 A 
E15 
LeS.5 8 
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224 58 
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E5 
44 6 reed bes 
Lae 
40.3 4.2 
E8 
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E9 
38 A 
gq=1.717cc/sec. 
RPM=250 
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A. 7 Ay 2 
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TABLE E-2 (continued) 
Xx10° sxi0"° x10" Sel” xeon Selo” seeLo- eo 
Inch Inch Inch Inch Inch Inch Inch Inch 
T3255 Pens E7 14 6.5 
17 12.9 6 6.3 E15 q=6.82cc/sec. 
E13 Tee addy Lee 5.3  RPM=1000 
8.5 16.1 E10 E17 El 
E14 1285 8.1 Re 34.7 30 6.2 
9 PSs7, 24 Gat. 2a @.5 18 8.6 
i3 15.2 31.5 4.0 E19 ies 8.6 
E15 Roll No. 149 3 2342-7  ROleeNeLOTEs 
a fy 13.5 E12 d=19.5x10 “in El 
E17 20 aie gs) E17 19.5 138 
16 14 E14 4 73 E2 
E18 aby ay TOUL OROLE No. 149 2.4 4.8 
13 25.27 25 9 E3 Beis HOw3 
15 pee der E15 6.5 no ee EA 
16.4 1258 § 29:25 fie La 10.8 94 16.8 
Roll No. 169 E16 E6 E9 
E5 8.5 6.7 18 se Sa ee ahe8 
11.6 16.6 Buz NCL=40x10 “in. 
E7 27 7123 @23.42cc/cees E1-E9 
3.9 np ee E18 RPM=500 2a #8 
NCL=40x1073in. 22 10.2 ‘Rol No. Dee 23.6 ze 
El NGL=20310-3in. 3.8 3 
6.5 4.06 g=1.71%ec/sec. Ell 18.5 8.3 
by 8.5  RPM=250 21.4 7.8 75 4.2 
E2 El EL2 32 5.9 
13,3 6:h 8692054 t6.8 3073 6.8 6 2.8 
24 Gud 23.5 8.8 EL oI pea, SAT 
E3 E3 22.5 7 Fie 742 
ea 4.667 12 8.5 E15 6 At 
20.7 Tote! 29 7.5 mete 5.6 10.7 5.2 
26.4 6.9 E6 E16 14 33 5.8 
31.6 Sa 29 16. “BBS 6.5 20.5 Se], 
38 3.04 Roll No. 169 E17 eG 4.9 
E5 E12 12.4 6,8 1805 5.6 
16.5 Sie oes Tilt) a6 a.9  2iea 8.9 
26 7.8 E14 E20 10 4.5 
E6 10.43 Tee lee Giese ae 
19:35 9.1 32.8 6.9 WNCL=3582%10-3in. 
E10 NCL=40x1073in. E13 Position No. 18 
15.5 $8 Bie 19.5 8.4 
27 6.3 26.5 6.9 Roll No. 148 Roll No. 170 
NCL=35.2x1073in. E13 E14 q=0.685cc/sec. 
E5 30 6.3 Weis 9.8 RPM=100 _, 
20.5 Bey E14 on 8.8  d=26.6x10 ~in. 
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TABLE E-2 (continued) 
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xx10° Sx10° xx10° Sx10? xx10° sxl0? xx10? Sx10° 
Inch Inch Inch fneh ineh £nen Inch Inch 
ELL q=1.717cc/sec. 23 12.9 533 12.8 
4.3 Loew RPM=250 Ez 
20.8 Ly a2 E14 Position No. 19 8.3 LY 
Baz 225.5 14.5 16.5 14.6 
10 IA AP E18 Row i Now 173 E14 
13 25 aL Gre) 20.8 q=0.685cc/sec. 5 11.9 
E13 Roll No. 170 RPM=100 3 AG pa oS 
13.5 2453 El d=2005%10°° in. NCL=27x1073in. 
se aay | 22 er. es: if wt ‘Beak E22 
E14 E2 1532 16.7 923 333 
v4 Lae 5a 6 14.6 E5 E14 
9,8 19.8 E3 3 OS L237 13,6 
L642 22.6 i2 iL? .'3 E7 NCL=50x1073in. 
E15 20 16.8 9.5 eke: ElL2 
8.8 1/23 2Y..2 LF L3z4 19.8 18 6.8 
18.8 23:28 E4 Luss LO«5 E14 
E17 iM) 20.5 E9 11.6 Dasa 
22.0 13 LS EOS 2eak (PS) 22.4 6.8 
E19 10 Mie rs: 
Na Ae | 225 q=3.42cc/sec. E10 q=6.82cc/sec. 
E20 RPM=500 L6<5 315-9 RPM=1000 
Dw LS Roll Noe L/7t NGL=27xL0NGin2 E3 
Sian J ro 2 El E9 Zoe sige 
ROLL NOG“ 150 133 4222 2.2 Le Ga) Sup 6.2 
E2 6.3 17-1 N€L=50x10/ San. E5 
BaD 203 12.6 24.8 225 rojas 20 38-9 
14 Zils iW 2360 d=2025x10 “ane 
La .D 2ka3 E4 q=1.717cc/sec. Bi 
19 L622 12 22.1 RPM=250 ~3 10.4 18.4 
E4 E5 d=20.5x10 ~in. E5 
Grek 2223 20.3 20.4 Ell 9 14.7 
be 19.4 L632 LO. J ig a ar sey 
16 20.3 gq=6.82cc/sec. EZ E6 
25 S23 RPM=1000 95 ieee" Loess Uh ore 
E7 BL2 18 te: Pee L236 16.8 
20 16.8 8 16.9 Ble 
23.2 6.03 EL3 17 12.5 Position No. 20 
E8 ] 16.6 E14 
5 LG.s 14 2576 Lb LV ARS: Refl No. 1524 
a ee 13.5 E15 q=0.685cc/sec. 
E10 EWA) 13.6 G=sin2cc/secr RPM=100 _, 
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i? .d 22.3 fers) LG.2. ‘Roll Noerl/2 E2 
15.5 28 Bil 93 LJ 
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TABLE E-2 (continued) 


Xx10° Sxl0° S10" Belo” xi0° Beto Sib” Sx10° 
Inch Inch Inch Inch Inch Inch Inch nen 
8.5 Vu staat 25.6 8.35 #1Lie3 5.8 53 9.4 
E5 42.6 4.2 45 S22 5452 9.1 
6.8 Te, 56.3 9.4 
16 10.8 q=1.717cc/sec. q=6.82cc/sec. 5 ia ecb 9.4 
E7 RPM=250 RPM=1000 602 9.3 
7 14.8 E1-E10 Baw 75.4 8.8 
Roll No. 174 17.5 726 30 9.8 85 7.6 
Ell 2427 Sp E12 94.2 aie ales: 
12 18 34 8.4 23 8.5 E5 
E12 24 Tse E13 2155 7.6 
eo 2392 1787 6.8 PAY. 10.6 58 9.4 
6.6 15.9 2.5 ae E15 63.2 9.4 
17.3 14.6 14 6.7 44.6 ve? 66.6 Oral. 
E14 35 8 E16 69.2 8.8 
19 10.4 Loe 5.4 9 S fe 100 6.6 
E15 Al 4.4 274 9.8 105 6.2 
6.8 16.7 19%8 6 ¥20 =3. 4:2 109.5 5.3 
9.8 L7e7 NCL=38.6x1074in. d=21.4x10 “in. E6 
14 16.7 El E12 33 8.3 
E16 1274 9.9 oo 6 19.6 36 8.6 
EL?3 19.4 25 9.4 56 9.4 
E17 E6 P6sitiom Now 21 90 ‘Loe 
EPs5 i a 25 G2 E7 
NCL=38.6x10-3in. d=21.4x10-3in. Roll No. 157 6 3.5 
jag ak EL q=0.685cc/sec. -5 seg 
20.3 i keer 7 17.4 RPM=100 _ 40 Sarg: 
E14 14 oom d=1.29x10 “in. 61.3 8.8 
324 6.3 E2 El 76.6 a7 
E16 B)F2 Pr.) 37 8.8 94 feral gs 
20 9.6 6 ‘ies hey) Mes 9.1 98 6.8 
NEE=49.2x1073in. E3 85 226 Te3 20 1,6 
E11-E20 8.4 18.4 89 Les E9 
29.3 8.35 E2 +e ih AS 6.3 
38.2 659  q=3t42cc/sec: 30.8 710 Goleerd pubes 
spl 3.34 RPM=500 98.3 9 96.5 7! 
6.6 3°76. Roll noe 175 60.35 8.4 E10 
A274 4.6 E2 89 ee 9.3 5.4 
18.4 oe 4 Fee? Aes Lge 34 74 ow 
33 8.8  Ne@b=49.2x107 in. E3 64 9 
172i 7.25 E5 cw Es: 8.8 Gi. 5 a. 3 
30 io 1255 6.8 55.8 9 90 iG 
30 8.6 35 See 0. Lo ue, oe fhe ea 
30.6 8.35 E7 E4 98.6 6.9 
15 7.3 6.4 4.8 aL erz ROLL NO. 170 
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TABLE E-2 (continued) 
3 3 3 3 
Xx10 Sx10 Xx10 SxlO”~ Xx10 Sx10 Xx10 Sx20 
Inch Inch Inch inch tIneh inch Inch Inch 
E16 Bil Sie ie ON E16 
29.8 eres) 8 6 49 19 30%. 2 TO s:6 
32.3 aca 10.9 6 68 19.2 
116 eae 37 9.5 Position No. 23 
Bey, 64.3 9.5 q=1.7l/cc/sec. 
or aes GSeo 10.1 RPM=250 Roll No. 204 
ins 6.4 69.5 9.4 E14 q=0.685cc/sec. 
86 1PrSE LOT <7 6.9 35.23 UR haa RPM=100 ~3 
E18 E12 71 aes NCL=53x10 “in. 
45 9.2 a, 9.1 ELS ek 
83.6 He'D 7804 Ceie (2a. 8 Le Bice Zia 
91 dD 74 L/ 42.7 oe) 
E19 q=6.82cc/sec. 98 16.2 2 
52 9.4 RPM=1000 Ey. 34.2 Ee | 
ADS 6 B22 lah 47.5 Ue ha 50 Gn03 
age 8.1 59.5 9.8 E18 E14 
E20 LO sD op Bad Les7 TP 5.9 
60.8 9.6 82.3 Bee Beleaen 52 
66 9 89 8.7 q=3.42cc/sec. 
92.4 Ls oO 8.9 RPM=500 q=1.717cc/sec. 
ye ae 100 io9 Boll Now L78 RPM=250 
Bz El E16 
g=1.717cc/sec. 61.8 9 46 17.9 9.4 18 
RPM=250 BS 9 65.5 La 2 48 ees: 
Rolt No, = 157 2225 B.D D4. 19.9 BL / 
Ell E2 11.4 Cline 
69 8.9 Position No. 22 37 14.6 E18 
70 8.9 41 oS aes | as: ik 
Ace G.3 Roll No. 179 3 19.5 23:10 
lon2 5.6 q=0.685cc/sec. 12.5 10.3 24 22.9 
121.9 Ses RPM=100 = 28 AERP} 318 i 2524 
E12 NCL=100.5x10 ~in.52.3 AINS gene’: 44°53 aly ey, 
GF eS oe5 Be az. Las 9 
43.5 08. 58 6.6 q=3.42cc/sec. 
50 9.4 79.4 19 q=6.82cc/sec. RPM=500 
B2n3S 9.8 EZ RPM=1000 ROdae (NOsacoil 
69.4 Deak. 2.00 455 Ell E2 
74 8.8 47 Zue3 “Oa 3.8 ibe gas y Sr 
87 8 64.5 L@.t 80.2.6 16.4 E5 
9145 fas 81.8 17.5 Biz 44 a ae 
E3 36 J53 4 E8 
q=3.42cc/sec. 28 16.6 69 1234 Sis0 19.4 
RPM=500 E4 Ba 
Roll No. 156 19 14.9) 20 12.9 q=6.82cc/sec. 
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TABLE E-2 (continued) 
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xx107 sx10? xx10° Sx10° xx10° Sig” eho? © 6shHOe 
Inch Inch Inch Inch inch Inch Inch Inch 
RPM=1000 E20 19 9 24.8 6.5 
E15 10 20.9 35 9 E3 
35.4 22 45 9 2.5 9.4 
q=3.42cc/sec. E9 8.3 25:38 
Position No. 24 RPM=500 10 Bet’ 22:58 138 
Roll No. 182 15 7.9 NCL=66x1073in. 
Roll No. 183 E17 LT oF 7.9 El 
q=0.685cc/sec. 12 22. 255 9 2h2 9.9 
RPM=100_, Roll No. 203 28.5 9 42 9.1 
d=22x10 ~in. NCL=59x1073in. E4 
El BL5 q=1.717cc/sec. SEeZ LOR 
8.1 19.8 36.6 6.8 RPM=250 51 6 
14 24.1 Ell E5 
E2 Position No. 25 19.5 10.1 15 7-26 
10.8 733 z 31.3 LO =1 
16 15.8 Roll No. 185 E12 q=6.82cc/sec. 
E4 q=0.685cc/sec. 23.5 9.4 RPM=1000 
9 Hose REN LO00 E13 Ell 
15.3 16.2 d=26.4x10 “in. 43 72S 8654 6 
E9 El 55 5 E12 
14 19.8 3.6 1474 29 sees ail g33 
NCL=59x1073in. 13.3 23.8 d=26.4e80 -in. E14 
El E2 Ell 16 pad 
26.5 7.2 BOS DOm2 28 202 E16 
E9 28 16.9 ‘Wiz 9658 19 8.7 
34.4 Tx2 E3 E12 E19 
ia 22 10 19 30.8 9.1 
q=1.717cc/sec. 16.3 23-2 LE 2029 36 9.6 
RPM=250_., E4 £35 21.6 d=26.4x1073in. 
d=22x10 ~in. gee 10.8 E13 E14 
E14 22 Gee ae 29.22 Tied 1993 
20 3.62 E5 4.3 23.6 9 20.4 
E15 U7 a2 27) Tq 25 LD 25 1774 
7 L327 E6 22 15.8 E15 
E16 8.3 19.8 4.1 16.5 
Ee 2A O21 19.4 q=3.42cc/sec. 6.8 2255 
5.4 13.4 NCL=66x1073in. RPM=500 9.5 LORS 
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